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PREFACE 


The  prediction  of  deepening  effects  on  sedimentation,  documented  in  this 
report,  was  performed  for  the  US  Army  Engineer  District,  Charleston,  and  in¬ 
cluded  use  of  techniques  recently  developed  under  the  Improvement  of  Opera¬ 
tions  and  Maintenance  Techniques  (IOMT)  research  program  sponsored  by  the 
Headquarters,  US  Army  Corps  of  Engineers,  under  IOMT  Work  Unit  No.  31765, 
"Fine-Grained  Shoaling  in  Navigation  Channels."  This  study  was  conducted  in 
the  Hydraulics  Laboratory  (HL)  of  the  US  Army  Engineer  Waterways  Experiment 
Station  (WES)  during  the  period  January  1986  to  September  1988  under  the  gen¬ 
eral  supervision  of  Messrs.  Frank  A.  Herrmann,  Jr.,  Chief,  HL;  Richard  A. 
Sager,  Assistant  Chief,  HL;  William  H.  McAnally,  Jr.,  Chief,  Estuaries 
Division;  and  George  M.  Fisackerly,  Chief,  Estuarine  Processes  Branch. 

The  study  was  conducted  and  this  report  prepared  by  Mr.  Allen  M.  Teeter, 
Estuarine  Processes  Branch.  Mr.  Walter  Pankow,  Estuarine  Processes  Branch, 
assisted  in  the  preparation  of  this  report.  Mr.  Howard  A.  Benson,  Estuarine 
Processes  Branch,  was  the  field  engineer  for  the  field  data  collection  which 
preceded  this  study.  Field  technicians  who  collected  data  included 
Messrs.  David  Crouse  (retired),  Joseph  W.  Parman,  James  T.  Hilbun,  Samuel  E. 
Varnell,  Billy  G.  Moore,  John  T.  Cartwright,  Douglas  M.  White,  and  John  S. 
Ashley,  all  with  the  Estuaries  Division.  Mrs.  Clara  Coleman,  Estuarine 
Processes  Branch,  reduced  the  data  to  computer  files.  Mrs.  Marsha  C.  Gay, 
Information  Technology  Laboratory,  WES,  edited  this  report. 

The  Charleston  District  contact  persons  were  Messrs.  Lincoln  Blake, 
Robert  Billue,  and  James  Joslin. 


Commander  and  Director  of  WES  during  preparation  of  this  report  was 


COL  Larry  B.  Fulton,  EN. 


Technical  Director  was  Dr.  Robert  W.  Whalin. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multiply 

By 

To  Obtain 

cubic  feet 

0.02831685 

cubic  metres 

cubic  yards 

0.7645549 

cubic  metres 

feet 

0.3048 

metres 

pounds 

0.000454 

metric  tons 

pounds  per  cubic  foot 

16.02 

grams  per  litre 

miles  (US  statute) 

1.6093 

kilometres 

square  feet 

0.09290341 

square  metres 
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SCHEMATIC  NUMERICAL  MODELING  OF  HARBOR  DEEPENING  EFFECTS 
ON  SEDIMENTATION,  CHARLESTON,  SOUTH  CAROLINA 


PART  I:  INTRODUCTION 

Background 

1.  The  US  Army  Engineer  Waterways  Experiment  Station  (USAEWES) 
Hydraulics  Laboratory  (HL)  has  studied  the  effects  of  the  Cooper  River  redi¬ 
version  on  sediment  flushing,  mixing,  and  circulation  in  Charleston  Harbor, 
Charleston,  SC  (Teeter  1989).  An  associated  study  task  was  to  reexamine  the 
possible  effects  of  ongoing  Charleston  Harbor  navigation  channel  deepening  on 
dredging  requirements  for  the  harbor.  Previously,  physical  model  hydrodynamic 
results  were  used  by  WES  to  qualitatively  assess  the  effect  of  harbor 
deepening  on  shoaling  (Benson  1976) . 


(  Purpose 

2.  The  purpose  of  this  study  was  to  predict  shoaling  for  the  new  proj¬ 
ect  channel  at  -AO  ft  mean  low  water  (mlw)  and  a  daily  average  inflow  of 
A, 500  cfs  for  the  purpose  of  providing  information  with  which  the  US  Army 
Engineer  District  (USAED) ,  Charleston,  planners  can  assess  disposal  site 
capacity  requirements. 


Scope 


3.  The  scope  of  this  study  included  the  schematic  application  of  a 
numerical  model  using  existing  hydraulic  and  sediment  information  to  gage  the 
effects  of  a  5-ft  channel  deepening  on  shoaling  rates  and  distribution  in  the 
main  portion  of  the  harbor  channel.  The  schematic  numerical  modeling  approach 
was  selected  because  no  simpler  (analytical  or  steady  state)  shoaling  predic¬ 
tor  was  known  which  could  account  for  both  channel  deepening  and  altered  in¬ 
flow  conditions.  A  two-dimensional  laterally  averaged  model  called  FIBS 
(Fine-Grained  Bed  Sediment)  was  employed.  A  limited  verification  of  the 
schematic  model  was  performed. 

A.  The  effects  of  rediversion  on  channel  shoaling  rates  were  also 


tested,  incidental  to  other  tests.  Salinity*  intrusion  and  stratification 
were  considered  by  this  study,  but  were  based  on  previous  physical  model  re¬ 
sults.  Modeling  included  the  effects  of  circulation,  and  covered  the  main 
harbor  channel  and  upstream  river  reach  to  the  Pinopolis  Dam  (Figure  1). 

5.  This  report  summarizes  the  methods  and  procedures  used  and  the 
results  obtained. 


*  Unusual  terms  are  listed  and  defined  in  Appendix  E. 
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PART  II:  STUDY  DESCRIPTION 


Prototype  Conditions 

6.  Charleston  Harbor  is  formed  by  the  junction  of  the  Ashley,  Wando, 
and  Cooper  Rivers,  and  is  a  major  South  Carolina  seaport  (Figure  1).  Average 
tide  range  is  5.3  ft  (1.61  m)*  at  the  Customs  House  at  Charleston,  located  at 
about  river  mile  9.  Conditions  in  Charleston  Harbor  have  recently  changed  as 
freshwater  inflow  was  reduced  by  about  71  percent  after  the  1985  rediversion 
of  the  Cooper  River.  The  effect  of  the  rediversion  on  shoaling  will  be  a 
gradual  decrease  in  harbor  shoaling  over  perhaps  a  decade.  Presently,  a  5-ft 
channel  deepening  is  underway  in  the  navigation  project,  possibly  altering 
shoaling  conditions.  Therefore,  Charleston  Harbor  is  in  transition  with  re¬ 
spect  to  shoaling,  and  the  only  shoaling  information  available  is  the  dredging 
records  for  pre-rediversion  harbor  conditions.  Because  of  the  recent  changes 
in  harbor  conditions  caused  by  the  rediversion  of  the  Cooper  River,  this  study 
addressed  the  effects  of  both  reduced  inflow  (rediversion)  and  increased 
channel  depths  on  harbor  shoaling. 

Background  Studies  and  Information 

7.  The  main  sources  of  information  for  this  study  were  dredging  records 
(USAED,  Charleston,  1966  and  I.  B.  Kyzer**),  previous  sediment  characteriza¬ 
tion  studies  (USAED,  Charleston,  1954  and  Teeter  1989),  physical  model  data 
(Benson  1976  and  USAEWES  1957) ,  and  sediment  budget  studies  (Patterson  1983 
and  Teeter  1989)  .  Two  physical  models  have  been  used  in  various  harbor 
studies  by  USAEWES,  including  two  studies  which  addressed  channel  deepening 
and  rediversion . 

8.  In  the  most  recent  physical  model  study,  USAEWES  tested  the  effects 
of  a  -40  ft  mlw  channel  on  tides,  currents,  and  salinities  (Benson  1976). 

Tests  were  performed  with  15,600-  and  3,500-cfs  inflows,  and  with  -35  and 
-40  ft  mlw  channels.  The  difference  between  mlw  and  National  Geodetic 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  found  on  page  4. 

**  Meeting  with  Mr.  I.  Braxton  Kyzer,  US  Army  Engineer  District,  Charleston, 
Engineering  Division,  in  April  1986. 
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Vertical  Datum  (NGVD)  is  1.98  ft.  The  difference  between  mlw  and  the  National 
Geodetic  Vertical  Datum  (NGVD),  which  is  now  the  standard  vertical  datum,  is 
1.98  ft.  However,  because  previous  navigation  projects  used  mlw  as  the  ver¬ 
tical  datum,  this  study  uses  it  also.  The  difference  between  mlw  and  mean  sea 
level  is  2.5  ft  at  the  Customs  House  tide  gage  at  Charleston. 

9.  The  physical  model  study  of  harbor  deepening  found  a  slight  increase 
in  tide  ranges,  very  slight  increases  in  currents,  and  a  slight  decrease  in 
salinity  stratifications  near  the  upper  end  of  the  harbor  project  associated 
with  the.  channel  deepening.  The  limit  of  the  100-ppm  saltwater  intrusion 
length  was  predicted  to  advance  about  1  mile  after  the  channel  deepening.  The 
study  concluded  that,  based  on  changes  in  hydrodynamic  conditions,  changes  in 
shoaling  conditions  would  be  minimal  with  a  tendency  for  peak  shoaling  to 
shift  toward  the  upstream  end  of  the  navigation  channel. 

10.  The  earliest  USAEWES  physical  model  studies  examined  the  effects  of 
the  rediversion  on  harbor  mixing  and  shoaling  (USAEWES  1957).  Tests  of  flows 
between  2,500  and  5,000  cfs  indicated  that  the  harbor  abruptly  became  well- 
mixed  and  sediment  flushing  improved.  Model  shoaling  tests  found  that  reduc¬ 
tion  in  freshwater  inflow  from  15,600  to  3,000  cfs  (with  a  corresponding 
decrease  in  sediment  inflow  of  80  percent)  reduced  shoaling  by  92  percent 
(exclusive  of  entrance  shoaling).  The  channel  deepening  from  -30  to  -35  ft 
mlw  which  occurred  in  the  1940's  was  found  to  have  had  a  minimal  effect  on 
harbor  shoaling. 

11.  The  US  Geological  Survey  (USGS)  also  studied  the  effect  of  rediver¬ 
sion  on  shoaling  (Patterson  1983).  Dredging  records,  hydrographic  surveys, 
maps,  charts,  hydrologic  data,  unpublished  files,  and  knowledgeable  individ¬ 
uals  were  sources  of  information  for  that  study.  Information  on  inflow  sus¬ 
pended  sediment  concentrations  was  compiled.  It  was  predicted  that  mainte¬ 
nance  dredging  would  be  reduced  by  40  to  75  percent  as  a  result  of  the 
rediversion. 

12.  Some  dredged  sediment  returns  to  the  estuary  due  to  overflows, 
leaks,  or  disposal  practices.  Patterson  (1983)  estimated  this  quantity  of 
sediment,  termed  runback,  to  be  22  percent  of  the  dredged  sediment  for  the 
1965-1984  period.  Average  annual  dredging  for  this  period  was  4.81  million 
cubic  yards  per  year  exclusive  of  the  entrance  channel. 

13.  The  Cooper  River  rediversion  was  carried  out  in  late  1985.  USAEWES 
performed  field  data  collection  and  analysis  on  current  and  salinity 
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conditions  before  rediversion  and  over  a  range  of  inflows  from  4,500  cfs  to  a 
project  base  level  of  3,000-cfs  daily  average.  A  small  amount  of  suspended 
material  data  was  also  collected  in  the  lower  harbor:  hourly  over  tidal 
cycles  and  composites  over  lunar  days.  Suspended  sediment  concentrations 
varied  between  50  and  350  mg/1  at  river  mile  7  (0.75  depth)  over  a  neap  to 
spring  tidal  sequence  and  4,900-cfs  inflow  (Figure  2).  A  storm  increased  sus¬ 
pended  sediment  concentrations  to  about  1,350  mg/1  for  2  days  at  that  loca¬ 
tion.  Analysis  of  suspended  sediment  data  from  two  surveys  at  river  mile  8 
indicated  that  tidal  pumping  was  the  dominant  flux  component  and  could  operate 
in  either  the  upstream  or  seaward  directions  (Teeter  1989). 

14.  A  study  of  rediversion  shoaling  reduction  between  3,000-  and 
4,500-cfs  daily  average  inflows  predicted  a  4  percent  difference,  from  74  to 
70  percent  reduction  respectively,  based  on  a  sediment  budget  approach  (Teeter 
1989).  That  analysis  is  presented  in  Appendix  A. 

Numerical  Modeling  Methods 

15.  The  numerical  sediment  model  used  in  this  study,  FIBS,  was  an  en¬ 
hancement  to  the  numerical  hydrodynamic  and  salinity  intrusion  model  system 
(SIMS  1)  developed  for  the  Cooper  River  rediversion  study  (Teeter  1989).  The 
FIBS  model  has  the  capability  to  predict  flow  and  transport  in  an  estuarine 
setting,  and  stresses  fluid  mud  bed  processes.  The  model  numerically  solved 
laterally  averaged  dynamic  differential  equations  for  flow  continuity,  conser¬ 
vation  of  horizontal  momentum,  and  conservation  of  mass  (salt  and  sediment) 
over  the  interval  between  the  water  surface  and  the  channel  bottom,  and  along 
the  length  of  the  channel.  The  horizontal  momentum  equation  included  nonlin¬ 
ear  advection  and  quadratic  bottom  friction  terms.  Appendix  B  describes 
details  of  the  flow  and  salt  transport  modeling. 

16.  The  basic  laterally  averaged  two-dimensional  flow  and  salt  transport 
equations  have  been  presented  previously  by  various  authors.  Such  systems  of 
equations  have  been  shown  to  reproduce  estuarine  circulation  driven  by  density 
currents  and  tidal  propagation.  In  FIBS,  additional  terms  were  included  in 
the  conservation  of  sediment  equation  to  account  for  sediment  dynamics.  Sedi¬ 
ment  dynamic  processes  for  settling,  deposition,  and  erosion  were  described  by 
a  set  of  algebraic  equations  which  depended  on  bed  shear  stress,  suspended 
sediment  concentration,  bed  sediment  concentration  (or  density)  and 
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sediment-dependent  coefficients  and  will  be  discussed  later. 

17.  Density  terms  were  included  in  the  momentum  equation,  coupling  the 
salt  transport  and  momentum  equations.  The  standard  assumptions  of  incom¬ 
pressibility,  hydrostatic  pressure,  and  negligible  vertical  momentum  were 
adapted,  allowing  vertical  current  velocities  to  be  computed  using  a  form  of 
the  continuity  equation. 

18.  Transport  equations  for  turbulent  kinetic  energy  with  production 
and  dissipation  terms  were  also  solved  in  the  model.  The  vertical  distribu¬ 
tion  of  eddy  coefficients  varied  according  to  an  analytic  expression.  For 
simplicity,  turbulence  closure  to  momentum  and  transport  equations  was  not 
implemented  for  this  study.  Vertical  distributions  of  eddy  diffusion  were 
used  in  the  calculations  of  effective  eddy  diffusivity  coefficients. 

19.  Vertical  eddy  dif fusivities  for  mass  and  momentum  were  functions  of 
input  coefficients,  depth,  and  Richardson  number.  Vertical  density  and  veloc¬ 
ity  gradients  were  used  to  dynamically  calculate  Richardson  numbers,  and  ad¬ 
just  vertical  eddy  dif fusivities  for  mass  and  momentum  according  to  an 
analytical  expression. 

20.  An  ocean  boundary  condition  was  formulated  for  FIBS  to  minimize 
reflection  of  long  waves  (surges)  originating  within  the  flow  domain.  The 
nonreflecting  model  ocean  boundary  permitted  the  dynamic  equilibrium  or  re¬ 
peating  tide  condition  to  be  attained  much  more  quickly  than  with  an  ocean 
water-level  boundary  specification  fixed  at  each  time  level.  Dynamic  equilib¬ 
rium  was  reached  in  about  two  tidal  cycles. 

21.  FIBS  solves  equations  over  finite  elements  using  a  method  of 
weighted  residuals  (orthogonal  collocation).  The  flow  domain  was  discretized 
as  a  series  of  elements  consisting  of  15  computational  points,  configured  as 
five  horizontal  rows  of  three  points  (or  nodes).  Each  node  consisted  of  a 
vertical  column  of  five  computational  points.  Breadth  at  each  node  varied 
vertically  as  a  quadratic  function  of  depth.  Model  equations  were  solved  im¬ 
plicitly  over  elements.  Time  integration  was  performed  using  an  explicit 
method  described  in  Appendix  B. 

22.  The  sediment  algorithm  for  FIBS  was  recently  developed  at  USAEWES  as 
part  of  a  research  work  unit  on  fine-grained  sediment.  A  number  of  associated 
algebraic  equations  were  solved  at  each  time-step  for 
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a.  Settling  velocity  W  *  of  sediments  as  a  function  of  suspension 
concentration  C  . 

b.  Depositional  probability  P  as  a  function  of  bed  shear  stress 

and  critical  shear  stress  for  deposition  t  ^  . 

c.  Critical  shear  stresses  for  particle  t  and  significant 
erosion  t  ^  as  functions  of  the  bed  density  Cg  . 

d.  Particle  and  significant  erosion  as  functions  of  critical  and 
bed  shear  stresses. 

Definitions  for  many  of  these  terms  and  details  on  sediment  transport  modeling 
are  presented  in  Appendix  C. 

23.  The  FIBS  model  accounts  for  changes  in  bed  concentration  or  density 
and  thickness  of  a  fluid  mud  layer  at  the  surface  of  the  bed.  The  bed  was 
defined  as  the  elevation  where  horizontal  motion  was  zero,  and  where  sediment 
concentration  was  greater  than  a  specified  value  (70  g/jt).  The  fluid  mud 
layer  was  defined  as  that  top  part  of  the  bed  where  sediment  concentrations 
were  less  than  the  fully  settled  condition  (300  g/ Z) .  Two  additional  bed- 
sediment  layers  were  included  in  the  model  below  the  fluid  mud  layer. 

24.  Formation  and  densif ication  of  the  fluid  mud  layer  were  computed 

using  differential  equations  for  continuity  and  conservation  of  mass,  and 

solved  over  the  variable  fluid  mud  vertical  domain  below  the  elevation  of  the 

bed.  A  numerical  scheme  similar  to  that  used  for  the  flow  domain  was  used  to 

solve  over  the  vertical,  five-computational-point  fluid  mud  layers  at  each 

node.  An  algebraic  expression  was  used  to  calculate  settling  rates  W,  as  a 

h 

function  of  bed  concentration  in  the  hindered-settling  range  of  concentrations 
C  .  Details  on  the  fluid  mud  model  are  given  in  Appendix  D. 


*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  defined 
in  the  Notation  (Appendix  F) . 
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PART  III:  TEST  CONDITIONS 


25.  Test  conditions  consisted  of  two  groups:  (a)  channel  and  inflow 
conditions  and  (b)  tide  and  sediment  conditions. 

Channel  and  Inflow  Conditions 


26.  Schematic  numerical  modeling  was  employed,  using  limited  spatial 
resolution  appropriate  for  system  definition  in  which  the  sediment  information 
was  limited.  The  numerical  mesh  covered  the  harbor  from  river  mile  1  to  river 
mile  25  at  1-mile  cross-section  intervals,  and  from  river  mile  25  to  Pinopolis 
(river  mile  57)  in  2-mile  cross-section  intervals  (Figure  1).  Figure  3  shows 
the  cross-sectional  areas  at  locations  along  the  estuary  and  river.  The  nu¬ 
merical  mesh  describing  estuarine  geometry  was  not  extensive  and  did  not 
include  branches  to  represent  the  East  Branch  Cooper  River  or  the  Back  Fiver 
Reservoir. 

27.  During  the  period  1965-1984,  the  project  channel  was  -35  ft  mlw, 
plus  -4  ft  of  advanced  maintenance  and  -2  ft  of  allowable  overdepth  dredging 
for  a  maximum  channel  depth  of  -41  ft  mlw.  The  Navy  channel  above  the  ordi¬ 
nance  reach  (river  miles  20-24)  had  a  maximum  depth  of  -37  ft  mlw  (-35  ft 
project  depth  plus  -2  ft  advanced  maintenance).  Since  1985,  overdepth  dredg¬ 
ing  was  reduced,  and  the  maximum  channel  depth  was  -39  ft  mlw.  All  model 
channel  depths  were  specified  as  2  ft  less  than  the  maximum  depths  to  repre¬ 
sent  average  channel  depths  during  a  dredging  cycle.  The  following  tabulation 
shows  channel  depths  for  the  three  harbor  conditions  modeled: 


Depths,  ft  mlw 


Condition 

Inflow 
cf  s 

Prototype 

Model 

Channel 

Proj  ect 

Maximum 

Harbor 

Navy 

Base 

15,600 

-35 

-41 

-39 

-37 

Base 

4,500 

-35 

-39 

-37 

-35 

Plan 

4,500 

-40 

-44 

-42 

-40 

Tide  and  Sediment  Conditions 


28.  Table  1  describes  seven  tidal  and  sediment  conditions  varied  for 
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test  prediction  and  sensitivity  analysis.  Table  1  also  defines  the  five  com¬ 
binations  of  tidal  and  sediment  conditions  used  for  test  predictions.  These 
are  discussed  in  Part  V. 

29.  Other  sediment  parameters  were  constant  for  all  tests.  For  the 
15,600-cfs  base  model  test,  initial  suspended  sediment  concentration  was 
250  ppm.  Suspended  sediment  concentrations  at  the  ocean  and  upstream  bound¬ 
aries  were  set  at  165  ppm  and  14  ppm,  respectively.  The  latter  is  the  average 
suspended  sediment  concentration  at  the  Pinopolis  inflow  to  the  Cooper  River, 
according  to  USGS  data. 

30.  For  4,500-cfs  base  and  plan  model  tests,  initial  suspended  sediment 
concentrations  were  75  ppm.  Suspended  sediment  concentrations  at  the  ocean 
and  upstream  boundaries  were  set  at  50  ppm  and  14  ppm,  respectively. 

31.  For  lack  of  site-specific  information,  known  sediment  characteris¬ 
tics  of  similar  estuarine  sediments  were  used  as  guidance  in  selecting  the 
necessary  sediment  parameters  for  this  study.  Preliminary  sensitivity  tests 
on  sediment  parameters  were  also  used  to  guide  sediment  parameter  selection. 

General  Test  Procedures 


32.  Model  test  sequences  were  used  to  represent  the  prototype  situation 
where  greater  than  average  concentrations  of  sediment  materials  are  in  suspen¬ 
sion  and  generally  erode  during  spring  tidal  conditions,  while  less  than  aver¬ 
age  concentrations  of  sediment  materials  are  in  suspension  and  generally  de¬ 
posit  during  neap  tidal  conditions.  Net  sedimentary  conditions  are  erosional 
between  neap  and  spring  tide  ranges,  and  depositional  between  spring  and  neap 
tide  ranges.  During  these  sequences,  suspended  sediments  are  added  upstream 
and  exchanged  with  the  ocean,  depending  on  sediment  flushing  conditions  in  the 
harbor.  Prototype  suspended  sediment  loads  thus  constantly  change  over  tides 
and  over  lunar  months,  and  the  balance  between  erosion  and  deposition  leads  to 
net  shoaling  or  erosion  over  time. 

33.  The  FIBS  numerical  model  was  applied  in  test  sequences  of  three 
steps  to  represent  prototype  suspended  sediment  variability. 

a.  The  first  computational  step  was  used  to  start  model  hydro¬ 
dynamics  and  gain  repeating  tidal  conditions  in  the  model.  An 
M2  lunar  tide  was  applied  at  the  ocean  boundary,  the  appropriate 
freshwater  inflow  was  introduced  upstream  at  Pinopolis,  and 
model  computations  were  carried  out  for  two  tidal  cycles. 


13 


b.  The  second  model  computational  step  was  an  erosional  phase  of 
one  tidal  cycle  produced  by  decreasing  fluid  mud  density  to 
allow  erosion. 

c.  The  third  and  final  computational  step  was  a  two-to-four  tidal 
cycle  equilibrium  phase.  Shoaling  rates  were  determined  at  the 
end  of  the  test  sequences. 
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PART  IV:  RESULTS 


Verification  Results 


34.  Model  tidal  propagation  was  verified  for  the  15,600-  and  4,500-cfs 
base  conditions  using  physical  model  results.  The  sediment  transport  model 
was  verified  to  observed  field  shoaling  rates  under  the  conditions  of 
15,600-cfs  inflow  and  -35  ft  mlw  project  channel.  The  average  annual  mainte¬ 
nance  dredging  records  from  1965-1984  were  used  for  the  shoaling  verification. 
Dredging  at  Shipyard  Creek  and  the  entrance  channel  was  not  included  in  the 
verification. 

Tidal  hydraulics 

35.  Tables  2-4  show  physical  and  numerical  model  tide  results  along  the 
main  channel  for  the  15,600-cfs  base,  4,500-cfs  base,  and  4,500-cfs  plan  con¬ 
ditions,  respectively.  Figure  4  shows  the  tidal  envelope  swept  out  by  instan¬ 
taneous  water-surface  profiles  for  physical  model  tests.  Figures  5-7  show  the 
numerical  and  physical  model  tidal  envelopes  for  the  15,600-cfs  base, 

4,500-cfs  base,  and  4,500-cfs  plan  conditions,  respectively.  Tidal  reproduc¬ 
tion  in  the  schematic  numerical  model  was  acceptable. 

Salinity 

36.  Physical  model  data  were  used  to  drive  salinity  conditions  in  the 
numerical  model.  Previous  Cooper  River  studies  adopted  chlorinity  to  repre¬ 
sent  salinity  because  the  concern  at  the  limit  of  intrusion  focused  on  chlo¬ 
rides.  Figures  8-10  show  salinity  distributions  along  the  harbor  channel  in 
the  numerical  model  for  the  15,600-cfs  base,  4,500-cfs  base,  and  4,500-cfs 
plan  conditions,  respectively.  Physical  model  salinity  results  were  used  to 
initialize  and  update  salinity  distributions,  forcing  numerical  model  and 
physical  model  results  to  be  very  similar.  It  was  not  within  the  scope  of 
this  study  to  predict  the  effects  of  deepening  on  salinity  distributions.  The 
physical  model  has  been  the  best  tool  applied  to  predict  salinity  effects  of 
channel  deepening  and  was  relied  upon  for  this  study. 

Shoaling 

37.  Shoaling  verification  results  are  shown  in  Table  5  and  Figure  11. 
They  show  the  average  prototype  dredging  rates  along  the  harbor  channel  for 
the  1965-1984  period,  and  the  values  predicted  by  the  model  for  15,600-cfs 
inflow  and  -35  ft  mlw  channel  condition.  The  agreement  between  the  overall 
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model  and  prototype  shoaling  magnitudes  was  excellent.  Upstream  of  river 
mile  11,  the  general  shape  of  the  verification  shoaling  distribution  was  simi¬ 
lar  to,  but  smoother  than,  the  prototype  distribution.  Between  river  miles  6 
and  10,  the  model  shoaling  distribution  was  spatially  out  of  phase  with  the 
prototype  distribution  probably  due  to  the  lack  of  branching  in  the  model. 
Although  the  comparison  of  shoaling  distributions  was  rather  poor,  it  was  con¬ 
sidered  sufficient  for  this  study.  Heavy  shoaling  areas  in  the  model  occurred 
where  maximum  near-bed  velocities  were  upstream,  and  near  the  point  of  zero 
tidal-ave raged  flow  (the  null  point).  Maximum  near-bed  velocities  and 
tidal-averaged  flows  are  discussed  in  the  next  section. 

Test  Results 


38.  Vertical  circulation  and  tidal  pumping  greatly  affect  harbor  shoal¬ 
ing.  Circulation  reflects  density  and  geometry  effects,  and  often  acts  to 
trap  sediments  in  estuaries.  Figures  12-14  show  tidal-averaged  circulation  in 
the  harbor  area  for  the  15,600-cfs  base,  4,500-cfs  base,  and  4,500-cfs  plan 
conditions,  respectively.  Near-bed  tidal-averaged  velocities  were  very  low  at 
river  miles  9-23  for  15,600-cfs  base,  and  at  river  miles  7-12  for  4,500-cfs 
base  and  plan  conditions. 

39.  Maximum  near-bed  velocities  reflect  density  and  tidal  effects,  and 
indicate  the  probable  direction  of  tidal  pumping  of  suspended  material. 

Figures  15-17  show  tidal  maximum  near-bed  velocities  squared  (proportional  to 
maximum  shear  stress)  and  the  direction  of  maximum  velocities  for  the 
15,600-cfs  base,  4,500-cfs  base,  and  4,500-cfs  plan  conditions,  respectively. 
Note  that  reversals  in  the  maximum  velocity  squared  values  reflect  relatively 
small  differences  in  longitudinal  values  and  in  ebb/flood  values.  Tidal- 
average  suspended  sediment  concentrations  are  shown  in  Figures  18-20  for  the 
15,600-cfs  base,  4,500-cfs  base,  and  4,500-cfs  plan  conditions,  respectively. 

40.  Table  6  presents  normalized  shoaling  results  for  all  three  channel 
and  inflow  test  conditions  and  for  the  tidal  and  sediment  conditions  defined 
in  Table  1.  Ratios  were  based  on  sediment  mass  of  predicted  shoals,  and  ex¬ 
press  the  plan  condition  results  normalized  by  both  base  conditions.  Results 
were  rounded  to  two  decimal  places.  Sediment  mass  gave  a  more  reliable  indi¬ 
cator  of  shoaling  than  shoal  volume  because  shoal  densities  in  the  numerical 
model  were  unverified  and  are  not  expected  to  vary  in  the  prototype.  Annual 


16 


shoal  deposits  in  the  prototype  would  be  thicker  layers  than  those  in  the 
numerical  model,  and  shoal  densities  would  be  expected  to  be  nearly  equal  for 
similar  shoaling  conditions.  Model  shoaling  predictions  are  for  the  main 
channel  of  the  harbor  and  do  not  include  the  seaward  entrance  channel  nor 
Shipyard  Creek  areas.  However,  these  areas  would  be  expected  to  experience 
similar  shoaling  patterns. 

41.  Longitudinal  shoaling  distributions  for  the  verification  tidal  and 
sediment  conditions  and  the  three  channel  and  inflow  conditions  are  shown  in 
Table  7.  The  rediversion  is  predicted  to  result  in  a  shoaling  distribution 
which  is  centered  farther  upstream.  The  center  of  the  distributions  changed 
from  river  mile  14  for  pre-rediversion  conditions  to  river  mile  21  for  base 
and  plan  conditions  at  4,500-cfs  inflow.  Thus  shoaling  reductions  were 
greater  in  the  lower  reaches  of  the  harbor. 

42.  The  numerical  model  predicted  erosion  for  river  miles  7-11.  This 
was  based  on  erosion  of  a  relatively  thin  fluid  mud  layer  and  should  not  be 
extrapolated  to  deeper  sediment  layers,  which  have  much  greater  hydraulic 
shear  strength.  However,  this  erosion  area  served  as  a  sediment  source  to  up¬ 
stream  areas  in  the  model,  and  would  physically  offset  some  shoaling.  Overall 
results  were  therefore  not  adjusted.  The  fluid  mud  modeling  is  described  in 
Appendix  D. 

43.  Table  8  shows  the  results  of  sensitivity  tests  of  select  tidal  and 
sediment  conditions.  As  shown,  results  were  particularly  sensitive  to  suspen¬ 
sion  and  bed  settling  velocities,  and  to  the  significant  erosion  rate  con¬ 
stant.  Model  shoaling  rate  decreased  by  about  one-third  between  tidal  cycles 
2  and  4  (Table  8,  last  column),  implying  that  longer,  spring-to-neap  tidal 
computations  might  improve  model  reliability. 
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PART  V:  DISCUSSION 


44.  Tests  for  4,500-cfs  base  and  plan  shoaling  were  made  using  a  number 
of  tidal  and  sediment  conditions,  as  defined  in  Table  I.  Apparently,  more 
than  one  combination  of  sediment  characteristics  would  yield  an  acceptable 
verification  to  observed  shoaling  (dredging)  rates.  In  addition,  the  model 
sensitivity  to  plan  tests  was  different  for  different  combinations  of  sediment 
characteristics.  Therefore,  since  there  were  uncertainties  in  specifying 
sediment  characteristics,  five  model  predictions  made  with  different  combina¬ 
tions  of  tidal  and  sediment  conditions  were  averaged.  This  procedure  was  used 
to  test  harbor  deepening  effects  and  improve  the  reliability  of  the  predic¬ 
tion.  Variability  in  the  model  prediction  for  a  given  channel  and  inflow  con¬ 
dition  is  an  indicator  of  the  overall  reliability  of  the  prediction. 

45.  Table  6  presents  the  model  test  results  expressed  as  a  ratio  of  the 
sediment  mass  deposited  in  each  test  over  the  sediment  mass  deposited  in  the 
15,600-cfs  base  test.  The  verification  test  of  tidal  and  sediment  conditions 
(a)  resulted  in  a  shoaling  rate  prediction  of  0.07,  or  7  percent,  for  the 
4,500-cfs  base  condition,  or  a  shoaling  reduction  of  93  percent.  While  this 
is  about  the  same  reduction  as  predicted  for  the  3,000-cfs  inflow  tests  by  the 
original  physical  model  studies  (USAEWES  1957) ,  it  is  much  higher  than  more 
recent  predictions  (Patterson  1983  and  Teeter  1989).  These  studies  are  dis¬ 
cussed  in  paragraphs  11  and  14,  and  more  thoroughly  in  Appendix  A.  Because  of 
this  difference,  condition  (a)  was  not  the  only  indicator  used  to  predict 
channel  deepening  effects  on  shoaling.  Four  additional  conditions  were  se¬ 
lected  for  testing  and  resulted  in  smaller  ranges  of  shoaling  reductions  than 
condition  (a).  The  five  tidal  and  sediment  conditions  (Table  1)  are  represen¬ 
tative  of  only  partial  segments  within  the  total  parameters,  but  cover  the 
expected  range  of  rediversion  shoaling  reduction.  Predictions  using  condi¬ 
tions  (c)  and  (d)  indicated  a  decrease  in  shoaling  after  harbor  deepening 
which,  intuitively,  does  not  appear  to  be  correct. 

46.  The  pre-rediversion  density  structure  was  used  in  condition  (e) ; 
therefore  this  test  generally  had  the  same  hydrodynamics  as  in  the  15,600-cfs 
base  condition.  There  were  other  parameters  used  in  condition  (e)  that  should 
have  reduced  the  shoaling.  The  parameters  of  lower  tide  range,  lower  suspen¬ 
sion  settling  coefficient,  and  lower  critical  erosion  stress  coefficient 
should  have  effectively  reduced  the  shoaling  according  to  the  sensitivity 
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results  (Table  8).  Therefore,  the  small  amount  of  shoaling  reduction  produced 
in  condition  (e)  is  attributed  to  the  stratified  harbor  conditions  imposed  in 
the  test.  However,  condition  (e)  did  approximate  the  effects  of  channel 
deepening  for  pre-rediversion  conditions,  although  the  prediction  of  these 
shoaling  conditions  was  not  the  focus  of  this  study. 

47.  When  conditions  (a)  through  (d)  are  averaged,  the  shoaling  rate 
ratio  is  about  0.32  for  the  4,500-cfs  base  and  plan  tests.  This  relates  to  an 
average  shoaling  reduction  of  68  percent,  which  is  reasonably  close  to  the  re¬ 
duction  predicted  by  the  sediment  budget  method  described  in  Appendix  A. 

48.  Since  the  purpose  of  this  study  was  to  provide  disposal  site  capac¬ 
ity  requirements,  a  conservative  approach  was  used  by  taking  the  larger  num¬ 
bers  of  the  predicted  range.  Another  reason  for  the  conservative  approach  was 
the  lack  of  site-specific  sediment  information.  The  maximum  shoaling  increase 
for  conditions  (a)  through  (d)  between  the  4,500-cfs  base  and  4,500-cfs  plan 
tests  was  0.04  relative  to  the  15,600-cfs  shoaling  rate,  and  is  another  con¬ 
servative  indicator  of  shoaling  increase  due  to  channel  deepening.  The  mean 
and  standard  deviation  for  the  4,500-cfs  base-to-plan  shoaling  differences 
relative  to  the  15,600-cfs  shoaling  rate  for  conditions  (a)  through  (d)  was 
0.00  and  0.04,  respectively.  The  mean  incremented  by  the  standard  deviation 
represents  another  measure  of  the  upper  predictive  range.  Therefore,  the  con¬ 
servative  model  prediction  for  increased  shoaling  due  to  channel  deepening  is 
0.04  times  the  shoaling  rate  of  the  15,600-cfs  base  condition  or  0.15  times 
the  shoaling  rate  of  the  4,500-cfs  base  condition,  and  is  summarized  in  the 
next  section. 
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PART  VI:  CONCLUSIONS 


49.  Harbor  deepening  to  -40  ft  mlw  was  predicted  to  increase  harbor 
dredging  requirements  by  about  15  percent  for  rediversion  conditions,  or 
210,000  cu  yd  annually  at  4,500  cfs.  The  rediversion  channel  deepening  effect 
was  the  main  focus  of  this  study,  and  the  increased  shoaling  will  amount  to 
about  4  percent  of  the  pre-rediversion  condition.  The  following  tabulation 
summarizes  the  combined  results  for  all  three  channel  and  inflow  conditions. 


Average  Inflow 
_ cfs _ 

15,600 

4,500 

3,000 


Charleston  Harbor 
Dredging  Requirement, 

10^  cu  yd  per  year,*  for 
Channel  Depth,  ft  mlw 

-35 

-40 

4.81** 

6.37* 

1.43** 

1.64* 

1.27** 

1.46* 

*  Excluding  entrance  channel  dredging. 

**  Prototype  average  for  1965-1984. 

t  Prediction  based  on  shoaling  rate  ratio  for 
tidal  and  sediment  condition  (e)  which  had 
stratification  and  density  currents  as  for 
15,600-cfs  base.  See  Tables  4  and  6,  and 
paragraph  46. 

tt  Prediction  based  on  sediment  budged  analysis. 
See  Appendix  A. 

t  Prediction  based  on  channel  deepening  effects 
from  schematic  numerical  model.  See  Table  6  and 
paragraph  48. 


50.  Entrance  channel  shoaling  was  not  predicted.  However,  no  increase 
in  shoaling  is  expected  as  a  result  of  the  harbor  channel  deepening  in  the 
entrance  channel  based  on  extrapolation  of  model  shoaling  distributions. 

51.  The  centroid  of  the  shoaling  distribution  is  expected  to  move  about 
7  miles  upstream  as  a  result  of  rediversion  of  the  Cooper  River.  Numerical 
model  results  indicated  that  the  harbor  channel  from  river  mile  7  to  11  will 
be  self-maintaining  after  redlversion  and  channel  deepening. 

52.  Numerical  modeling  of  fine-grained  sediment  transport  in  estuaries 
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is  not  a  simple  matter.  The  model  did  not  respond  In  a  simple  way  to  changes 
in  channel  depth.  Results  depended  in  a  complex  way  on  other  tidal  and  sedi¬ 
ment  conditions  imposed.  Therefore,  numerical  hydrodynamic  and  sediment 
models  do  not  lend  themselves  to  simplified  application  to  these  types  of 
problems.  In  this  study,  sensitivity  testing  was  used  to  examine  a  range  of 
possible  effects,  and  a  value  at  the  conservative  end  of  the  range  selected  to 
characterize  the  response. 
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Table  1 

Tidal  and  Sediment  Conditions 


Tidal 

and  Sediment  Condition  Designation 

Tide  and  Sediment  Condition 

(a)* 

(b) 

(c) 

(d) 

(e) 

Tide  Range  at  Ocean,  m 

1.70 

1.70 

1.70 

1.70 

1.60 

Salinity 

Distribution 

Normal 

Normal 

Normal 

Normal 

15,600-cfs 

base 

Suspension  Settling 
Coefficent  A1 

(W  =  Al  C4/3) 
s 

1 . IE- 7 

1.4E-7 

1.6E-7 

1.8E-7 

1.0E-7 

Critical  Stress  for 
Deposition 

Tcd  »  N/sq  m 

0.06 

0.08 

0.08 

0.08 

0.08 

Fluid  Mud  Settling 

Coefficient  W1  ,  cm/sec 

(W  =  W1  •  f(C  )) 
n  s 

IE-4 

IE-4 

5E-5 

IE-5 

IE-4 

Critical  Erosion  Stress 
Coefficient  Dhl  , 

(t  ,  -  Dhl  •  C  5/2) 
ch  s 

90 

90 

90 

90 

40 

No.  of  Depositional 

Tidal  Cycles  (step  3) 

2 

4 

2 

2 

2 

*  Verification  condition. 

Table  2 


Tidal  Verification  for  15,600-cfs  Base  Condition 


River 

Model  Tide 

Range,  m 

Model  Average 

Phase,  hr 

Mile 

Node 

Physical 

Numerical 

Physical 

Numerical 

1 

1 

1.52 

1.73 

0 

0 

9 

11 

1.80 

1.81 

0.4 

0.3 

19 

21 

1.86 

1.87 

1.2 

1.0 

25 

27 

1.80 

1.91 

1.5 

1.3 

33 

31 

1.68 

1.67 

2.3 

2.0 

41 

35 

1.13 

1.64 

3.2 

2.5 

45 

37 

0.64 

1.17 

4.1 

3.3 

51 

40 

0.55 

1.12 

5.0 

3.7 

Table  3 

Tidal  Verification  for  4,500-cfs*  Base  Condition 


River 

Model 

Tide  Range,  m 

Model  Average 

Phase,  hr 

Mile 

Node 

Physical 

Numerical 

Physical 

Numerical 

1 

1 

1.68 

1.71 

0 

0 

9 

11 

1.83 

1.71 

0.3 

0.5 

19 

21 

1.80 

1.64 

1.1 

1.5 

25 

27 

1.74 

1.66 

1.4 

1.5 

33 

31 

1.68 

1.45 

2.1 

2.5 

41 

35 

1.55 

1.38 

3.0 

3.0 

45 

37 

0.94 

1.12 

3.9 

3.5 

51 

40 

0.94 

1.00 

4.1 

4.0 

*  Physical  model  had 

3,500-cfs 

inflow. 

Tidal 

Comparison 

Table  4 

for  4,500-cfs*  Plan 

Condition 

River 

Model  Tide  Range,  m 

Model  Average 

Phase,  hr 

Mile 

Node 

Physical 

Numerical 

Physical 

Numerical 

1 

1 

1.65 

1.72 

0.0 

0.0 

9 

11 

1.86 

1.76 

0.3 

0.3 

19 

21 

1.86 

1.79 

0.9 

1.0 

25 

27 

1.83 

1.81 

1.3 

1.5 

33 

31 

1.68 

1.55 

2.0 

2.2 

41 

35 

1.55 

1.51 

3.0 

2.7 

45 

37 

0.91 

1.19 

3.7 

3.5 

51 

40 

0.98 

1.08 

4.1 

4.0 

*  Physical  model  had  3,500-cfs  inflow 


Table  5 

Shoaling  Verification  for  15,600-cfs  Base  Condition 


River  Mile 

After 

Prototype  Dredged 

Volumes 

Verification 
Dredged  Volumes* 
cu  m/cycle 

Shoal/ 

Reach 

River  Mile 

Node 

1,000  cu  yd/year  < 

cu  m/cycle 

6 

8 

315 

319 

901 

Anchorage 

7 

9 

452 

458 

479 

Anchorage 

8 

10 

0 

0 

223 

— 

9 

11 

713 

723 

10 

Tidewater , 

6A,6C 

10 

12 

620 

628 

15 

6A,  6C 

11 

13 

134 

136 

223 

6A.6B.6 

12 

14 

58 

59 

261 

6B,6 

13 

15 

360** 

365** 

233 

5A 

14 

16 

46 

47 

209 

5A 

15 

17 

0 

0 

189 

— 

16 

18 

170 

172 

177 

4 

17 

19 

9 

9 

173 

3 

18 

20 

48 

49 

220 

1&2 ,3 

19 

21 

200 

203 

301 

1&2 

20 

22 

120 

122 

370 

1&2,  NAD 

21 

23 

400 

405 

303 

NAD 

22 

24 

378 

383 

83 

NAD 

23 

25 

10 

20 

5 

NAD 

TOTAL 

4,145 

4,088 

4,374t 

*  Dredged  volumes  =  1.28  times  shoaling  volumes,  assuming  22  percent 
runback. 

**  Excluding  Shipyard  River  which  amounts  to  780,000  cu  yd  annually, 
t  On  average,  Cs  =  0.334  g/cu  cm. 


Table  6 

Shoaling  Test  Results 


Shoaling  Rate 

Ratio*  for  Channel 

and  Inflow 

Condition 

Tidal  and** 

Relative 

to  15,600-cfs  Base 

Relative 

to 

4,500-cfs 

Base 

Sediment 

Condition 

15,600-cfs 

Base 

4,500-cfs 

Base 

4,500-cfs 

Plan 

4,500-cfs 

Plan 

(a) 

1.0 

0.07 

0.09 

1.33 

(b) 

0.17 

0.21 

1.22 

(c) 

0.A2 

0.41 

0.97 

(d) 

0.62 

0.57 

0.91 

(e) 

0.70 

0.93 

1.32 

Averages : 

0.40 

0.44 

1.15 

Standard  Deviations: 

0.28 

0.33 

0.20 

*  Ratio 


mass  of  sediment  deposited 


** 


mass  of  sediment  deposited  in  base  test  * 

See  paragraph  AO  in  text  and  Table  1  for  explanation. 


Table  7 

Model  Shoaling  Distributions 


Model  Shoaling  Distributions  for 
Verification  Conditions 

Mile  After _  _ Metric  Tons  per  Tidal  Cycle 


River 

Mile 

Node 

15 ,600-cf s 

Base 

4,500-cfs 

Base 

4 , 500-cf s 
Plan 

6 

8 

172 

5 

8 

7 

9 

109 

-5 

-1 

8 

10 

68 

-10 

-6 

9 

11 

35 

-15 

-9 

10 

12 

30 

-14 

-10 

11 

13 

55 

-4 

-4 

12 

14 

60 

0 

0 

13 

15 

56 

4 

3 

14 

16 

53 

7 

7 

15 

17 

50 

8 

8 

16 

18 

48 

6 

6 

17 

19 

47 

7 

8 

18 

20 

57 

11 

12 

19 

21 

71 

14 

13 

20 

22 

83 

15 

14 

21 

23 

79 

17 

16 

22 

24 

49 

20 

21 

23 

25 

21 

19 

28 

TOTAL 


1,142 


84 


122 


Channel  and  Inflow  Condition 


Shoaling  rate  924  &  361  924  &  3»947  924  &  863  480  &  8,447  305  &  193 

m.  tons  per  cycle 
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Figure  3.  Cross-sectional  areas  along  the  Cooper 
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Figure  5.  Water-surface  profiles  at  half-hourly  intervals  for  1 5,600-cfs  base  test 
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Figure  6.  Water-surface  profiles  at  half-hourly  intervals  for  4,500-cfs  base  test 
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Figure  7.  Water-surface  profiles  at  half-hourly  intervals  for  4,500-cfs  plan  test 
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Figure  II.  Verification  shoaling  distribution  for 
15,600-cfs  base  condition 
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Figure  13.  Profile  of  vertical  circulation  vectors  along  the  harbor  channel  for  4,500-cfs  base  test 


Figure  15.  Maximum  near-bed  velocities  squared  and  directions  along  the 

harbor  channel  for  15,600-cfs  base  test 


Maximum  near-bed  velocities  squared  and  directions  along  the  harbor  channel 
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Figure  19.  Profile  of  suspended  sediment  values  along  the  harbor  channel  for  4,500-cfs  base  test 
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APPENDIX  A:  SHOALING  RATES  AT  3,000-  TO  4,500-CFS  INFLOW 


Purpose 


1.  The  purpose  of  this  appendix  is  to  present  an  analysis  of  the  dif¬ 
ferences  in  Charleston  Harbor  shoaling  for  the  range  of  inflows  from  3,000  to 
4,500  cfs. 


Brief  Review  of  Previous  Shoaling  Analyses 

2.  Study  of  the  causes  of  the  Charleston  Harbor  shoaling  began  shortly 
after  the  diversion  of  the  Santee  River's  flow  into  the  Cooper  River  in  1942. 

A  model  study  at  the  US  Army  Engineer  Waterways  Experiment  Station  (USAEWES) 
was  authorized  in  1947,  and  complementary  studies  were  carried  at  the  US  Army 
Engineer  District  (USAED) ,  Charleston,  in  the  late  1940's.  By  the  1950's, 
experience  with  the  diversion  and  investigations  into  the  harbor  shoaling  had 
identified  the  diversion  as  the  major  cause  of  the  shoaling  problem. 

3.  Prototype  investigations  by  USAED,  Charleston  (1954),*  showed  the 
suspended  material  in  the  Cooper  River  to  be  identical  to  those  in  the  shoal 
material.  Kaolinite  was  identified  as  the  primary  mineral  in  the  fines.  Set¬ 
tling  tests  showed  that  75  percent  of  the  particulates  settled  at  less  than 
0.001  fps.  Shoal  densities  and  grain  size  distributions  were  measured.  By 
all  indications,  the  shoal  material  was  characteristic  of  fine-grained,  co¬ 
hesive  sediments.  During  freshet  conditions,  highly  turbid  waters  were  ob¬ 
served  to  pass  from  the  Lake  Moultrie  Reservoir  into  the  Cooper  River  and  fill 
the  harbor  with  a  reddish  hue.  Thin  laminations  of  s'ick  clayey  materials, 
notable  by  their  distinctive  color  and  texture,  formed  on  tidal  flats  and 
beaches  in  several  areas. 

4.  This  study  concluded  that  the  rivers  were  the  largest  initial 
source,  with  settling  occurring  in  areas  of  relative  stagnation,  trapping 
material  in  the  flood-dominated  bottom  levels  of  the  estuary.  Material  near 
trie  bed  can  move,  dependent  on  its  density  and  viscosity,  until  it  has  the 
opportunity  to  harden  sufficiently. 


*  References  cited  in  this  appendix  can  be  found  at  the  end  of  the  main  body 
of  the  report. 
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5.  The  earliest  model  study  (USAEWES  1957)  concluded  that  more  than 
99  percent  of  the  shoaling  increase  was  due  to  the  diversion,  brought  on  by 
the  following  factors: 

a.  A  density  flow  superimposed  on  the  tidal  flow  that  produced 
strong  flood-dominated  flows  near  the  bed,  preventing  the 
estuary  from  disgorging  its  load  to  the  sea. 

b.  Increased  colloids  and  dissolved  material  available  to  shoal 
the  harbor,  both  from  suspended  load  in  the  river  and  from 
erosion  of  the  upper  channel  (assumed  to  have  equal 
magnitudes) . 

Channel  deepening  from  -30  to  -35  ft  mean  low  water  (mlw)  was  found  to  have 
caused  only  a  minor  increase  in  harbor  shoaling.  Rediversion  model  tests 
performed  at  2,500-  and  5,000-cfs  inflow  implied  that  about  3,000  cfs  was  the 
maximum  tolerable  to  harbor  stratification,  but  the  report  cautioned  that  no 
single  flow  was  best  for  the  entire  harbor. 

6.  The  latest  study  of  shoaling  conditions  in  Charleston  Harbor  and  the 
effect  of  rediversion  on  shoaling  was  performed  by  the  US  Geological  Survey 
(Patterson  1983)  .  An  attempt  was  made  to  quantify  sediment  sources  for  the 
system,  and  balance  them  against  amounts  of  sediment  removed  by  dredging  and 
storage  in  deposits.  Patterson  gathered  existing  information  to  estimate 
rates  of  sediment  inflow,  removal,  and  accumulation.  Dredging  records,  hydro- 
graphic  surveys,  maps,  charts,  hydrologic  data,  unpublished  files,  and  knowl¬ 
edgeable  individuals  were  sources  of  information  for  this  study. 

7.  Patterson  divided  the  data  into  approximately  20-year  periods  to 
identify  trends.  The  sediment  sources  identified  (and  mean  annual  values  for 
1966  through  1982)  included  the  following: 

a.  Pinopolis  discharge  (0.8  million  cubic  yards). 

b.  Cooper  River  scour  (0.25  million  cubic  yards). 

£.  Background  sources  including  diatom  plankton,  marsh  vegetation, 
urban  storm  runoff,  wastewater,  and  shoreline  erosion  (0.2, 

0.6,  0.15,  0.02,  and  0.3  million  cubic  yards,  respectively). 

d.  Unknown  (ocean  and  unspecified)  sources  (3.4  million  cubic 
yards) . 

The  unknown  source  magnitude  was  estimated  by  subtracting  the  known  inputs 
from  the  total  of  the  amounts  removed  and  accumulated.  This  study  was  able  to 
account  for  only  less  than  half  of  the  shoaling  by  known  sediment  sources. 

The  Pinopolis  discharge  of  suspended  sediments  (which  was  estimated  using  a 
number  of  different  methods)  accounted  for  only  about  15  percent  of  shoaling 
volumes . 
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8.  By  the  early  1960’s  Charleston  Harbor's  shoaling  rate  had 
stabilized  with  respect  to  the  diversion.  Dredging  became  more  effective  in 
the  early  1950' s,  and  river  channel  erosion  greatly  diminished. 

9.  Sediment  sources  were  projected  for  3,000-cfs  inflow  to  predict  the 
effect  of  rediversion  on  shoaling  rates.  A  shoaling  reduction  of  40  to 
75  percent  was  predicted  for  the  Cooper  River  rediversion,  based  largely  on 
the  unknown  component . 

Shoaling  Processes  in  Charleston  Harbor 

10.  The  1942  diversion  resulted  in  an  increase  in  freshwater  inflow 
from  about  600  to  15,600  cfs,  and  caused  about  a  36-fold  increase  in  inner 
harbor  dredging  (from  120,000  cu  yd  per  year  in  1953  exclusive  of  bar  and 
jetty  channel)  as  well  as  substantial  increases  in  other  areas  (USAEWES  1957). 
Three  hydrodynamic  sediment  traps  were  created  by  the  diversion,  and  were 
largely  responsible  for  increased  retention  of  shoaling  material  and  buildup 
of  unconsolidated  mud  throughout  the  estuary: 

a.  Vertical  density  stratification  increased  drastically  and 
trapped  sediments  near  the  bed. 

b.  Net  tidal-averaged  circulation  patterns  changed  and  trapped 
near-bed  suspended  sediments  in  developed  areas  of  the  estuary. 

c.  Once  concentrated  and  deposited,  sediments  were  trapped  in 
unconsolidated  mud  and  isolated  to  a  large  extent  from 
transport  by  turbulent  tidal  flows. 

The  increase  in  project  channel  depth  that  occurred  at  about  the  same  time  as 
the  diversion  did  not  have  an  important  effect  on  harbor  shoaling. 

11.  The  major  effects  of  rediversion  on  shoaling  for  both  the  4,500- 
and  3,000-cfs  weekly  average  flows  were  as  follows: 

a.  Reduce  vertical  density  stratification,  thus  improving  vertical 
mixing,  preventing  sediments  from  being  trapped  near  the  bed, 
and  improving  sediment  flushing  for  the  harbor. 

b.  Move  the  null  area  of  vertical  circulation  upstream,  thus  re¬ 
ducing  suspended  sediment  accumulation  and  unconsolidated  mud 
formation  in  project  and  facility  areas. 

c.  Reduce  sediment  and  nutrient  loadings  to  the  harbor. 

Vertical  density  stratification,  which  is  created  by  vertical  salinity  strati¬ 
fication,  damps  vertical  mixing.  The  null  area  of  vertical  circulation  is 
where  near-bed  net  (tidal-averaged)  velocities  are  neither  landward  nor 


A3 


seaward,  and  thus  is  an  area  of  converging  net  bottom  flow.  Organic  materials 
contribute  to  shoaling  directly,  and  foster  the  coagulation  of  inorganic 
sediments. 

12.  The  largest  reservoir  of  potential  channel  shoaling  material  is  now 
the  vast  (20-  to  30-million-cubic-yard)  blanket  of  unconsolidated  mud  that 
covers  the  floor  of  the  estuary.  This  material  has  densities  between  1.22  and 
about  1.05  g/cu  cm,  and  consistencies  between  that  of  mayonnaise  and  pea  soup. 
Unconsolidated  mud  has  been  observed  to  move  within  estuaries.  It  can  move 
longitudinally  landward  or  seaward  in  response  to  changing  tidal  and  fresh¬ 
water  inflow  conditions,  or  laterally  due  to  channel  slopes  or  special  flow 
conditions.  These  sediments  are  not  generally  moved  with  the  net  estuarine 
circulation  as  are  suspended  sediments.  Unconsolidated  muds  slump  or  move 
only  with  stronger  tidal  flows  near  the  bed,  and  tend  to  accumulate  in  deeper 
areas  of  relative  stagnation. 


Prediction  Method 


13.  Prediction  of  the  difference  in  shoaling  between  3,000  and 

4,500  cf 8  was  made  using  a  method  similar  to  that  used  by  Patterson  (1983).  A 
sediment  budget  was  constructed  for  Charleston  Harbor  that  identified  various 
sediment  source  components.  Data  on  average  annual  sediment  sources  from  1966 
through  1982  were  used.  The  effect  of  rediversion  on  each  component  was  esti¬ 
mated  to  make  shoaling  predictions.  Overall  shoaling  for  the  Charleston 
Harbor  was  considered. 

14.  The  latest  Charleston  District  estimate  of  the  1965  through  1984 
average  annual  gross  dredging  for  Charleston  Harbor  (6.19  million  cubic  yards 
per  year)  was  used  in  this  analysis,  and  was  somewhat  lower  than  Patterson's 
value  of  7.6  million  cubic  yards  per  year.  The  Charleston  District  value 
Includes  the  Naval  Ammunition  Depot  (NAD)  channel,  shoals  1-6,  Customs  House 
and  tidewater  reaches,  Shipyard  River,  anchorage,  and  entrance  channel 
dredging.  Annual  dredging  rates  and  locations  of  these  major  shoals  are  given 
in  Table  Al.  It  does  not  include  Navy  and  other  slips,  which  amounted  to 
3.13  million  cubic  yards  per  year  average  for  the  period  1953  through  1963. 

The  pier-slip  contribution  to  1965  through  1984  dredging  is  not  known,  but  is 
probably  lower  than  the  older  average.  The  estimated  difference  in  shoaling 
rates  between  3,000  and  4,500  cfs  was  not  affected  by  this  omission.  Runback, 
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the  difference  between  gross  dredging  volumes  and  permanent  removal,  was 
assumed  to  be  22  percent. 

15.  Plant  contributions  to  shoaling  by  marshes  and  diatom  plankton  were 
treated  separately  in  this  analysis,  as  they  are  considered  to  depend  on 
inflow.  The  high  level  of  productivity  in  estuarine  and  coastal  waters  has 
been  attributed  in  part  to  enrichment  by  nutrients  carried  by  river  waters 
(Parsons  and  Takahashi  1973) .  Additional  estuarine  biological  productivity 
enhancement  comes  from  the  mixing  of  fresh  water  (in  which  phosphorus  limits 
plant  growth)  with  ocean  water  (in  which  nitrogen  usually  limits  plant 
growth),  and  from  the  entrainment  of  deeper,  nutrient-rich  coastal  waters  by 
estuarine  flows.  Plant  production  contributions  were  assumed  to  come  from 
dissolved  and  particulate  nutrients,  largely  nitrogenous  materials,  carried  by 
the  inflow. 

Expected  Shoaling 

16.  The  difference  in  direct  sediment  inflow  and  plant  production  be¬ 
tween  4,500-  and  3,000-cfs  weekly  average  flows  will  amount  to  about 
160,000  cu  yd  of  shoaling  material  annually  (Table  A2) .  Sediment  inflow  and 
plant  production  contributions  to  shoaling  are  expected  to  be  proportional  to 
Pinopolis  inflow. 

17.  The  unknown  sediment  source  referred  to  in  Table  A2  could  be  made 
up  largely  from  sediments  of  ocean  origin.  Reduction  of  the  unknown  source 
was  related  to  the  improved  sediment  flushing  efficiency  of  the  harbor,  and 
therefore  inversely  proportional  to  the  vertical  density  stratification  ob¬ 
served  during  flow  testing  surveys.  Scour  in  the  Cooper  River  is  expected  to 
be  eliminated  for  both  3,000-  and  4,500-cfs  flows. 

18.  The  overall  shoaling  reduction  predicted  in  Table  A2  for  3,000  cfs 
(74  percent)  is  slightly  greater  than  the  Charleston  District  overall  1966 
estimate  (71  percent),  and  slightly  less  than  the  upper  limit  of  Patterson's 
(1983)  predicted  range  (40-75  percent). 

19.  The  overall  difference  in  dredged  volumes  between  4,500-  and 
3,000-cfs  weekly  average  flow  will  most  likely  be  about  200,000  cu  yd 
annually. 
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Entrance  Channel  Shoaling 


20.  The  overall  shoaling  estimates  presented  in  the  preceding  para¬ 
graphs  included  the  entrance  channel.  Entrance  channel  shoaling  will  be  con¬ 
sidered  separately  in  this  section  because  of  dredging  cost  concerns  for  this 
area  and  because  of  the  paradoxical  nature  of  shoaling  seaward  from  a  harbor 
many  experts  have  declared  to  be  an  efficient  sediment  trap. 

21.  Early  studies  described  Charleston  Harbor  entrance  channel 
shoaling  material  as  coarse-grained.  Patterson  (1983)  compiled  prediversion 
dredging  volumes  for  the  entrance  channel.  Long-term  averages  were  less  than 
267,000  cu  yd  per  year.  Recently  the  Charleston  District  has  identified 
entrance  channel  shoaling  as  fine-grained.  Average  entrance  channel  dredging 
was  1.24  million  cubic  yards  per  year  for  1965  through  1984. 

22.  Rough  calculations  of  deposition  rates  from  sediment  suspensions 
were  made  for  the  entrance  channel  using  reasonable  values  for  settling  veloc¬ 
ity  (0.01  cm/sec),  near-bed  suspension  concentration  (200  mg 1 1) ,  and  the  fre¬ 
quency  of  deposition  time.  The  frequency  of  deposition  time  was  estimated 
using  an  assumed  critical  shear  stress  for  deposition  (0.05  Pa)  and  compiled 
coastal  currents  for  the  area  (0.3  to  0.8  knots)  as  8  percent  (about  30  days 
per  year).  To  balance  the  observed  shoaling  mass  (specific  weight  of  shoal 
material  times  shoal  volume  divided  by  shoal  area)  with  calculated  deposition 
from  suspension  required  unreasonable  values  for  depositional  frequency 

(267  days  per  year)  or  for  near-bed  concentration  (1,780  mg / £) .  Therefore,  it 
is  difficult  to  account  for  shoaling  in  the  entrance  channel  by  settling  from 
suspension,  even  when  the  possibility  of  reerosion  of  deposited  sediments  by 
storm  action  is  totally  ignored.  Therefore,  even  when  the  possibility  of  re¬ 
erosion  of  deposited  sediments  by  storm  action  was  ignored,  only  a  small  part 
of  the  shoaling  in  the  entrance  channel  was  attributed  to  settling  from 
suspension. 

23.  Entrance  channel  shoaling  increases  have  probably  been  caused  by 
near-bed  movements  of  unconsolidated  mud.  Ebbing  tidal  flows  transported 
sediments  out  of  the  estuary,  and  they  became  stranded  in  the  outer  entrance 
channel  where  tidal  flood  flows  were  insufficient  to  return  them. 

24.  It  is  reasonable  to  assume  that  there  will  be  a  considerable  reduc¬ 
tion  in  entrance  channel  shoaling  after  rediversion  and  a  stabilization 
period.  This  was  also  the  opinion  of  the  Committee  on  Tidal  Hydraulics 
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(1966).  The  shoaling  reduction  in  this  area  will  be  of  the  same  order  of  mag¬ 
nitude  as  the  predicted  overall  reduction.  The  unconsolidated  mud  shoaling 
source  will  diminish  over  the  next  decade — some  of  it  dredged  from  channel 
sites,  some  flushed  seaward  from  the  harbor,  and  some  hardening  in  place. 
Suspended  sediment  flushed  seaward  is  not  expected  to  deposit  rapidly  enough 
to  increase  entrance  channel  shoaling. 
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Table  A1 


Annual  Dredging  Rates  and 
Locations  of  Major  Inner  Harbor  Shoals,  1965-1984* 


Shoal/Reach 

Annual  Gross 

Yardage 

1,000  cu  yd 

River  Mile 

Anchorage 

767 

6.5-  7.7 

Tidewater 

563 

9.1-  9.7 

6A 

547 

9.7-11.1 

6B 

37 

11.4-12.2 

6C 

263 

9.9-10.7 

6 

115 

11.6-12.3 

Shipyard  River 

780 

13.0-13.7 

5A 

406 

13.2-14.1 

4 

170 

16.2-17.0 

3 

37 

17.7-18.9 

1  and  2 

300 

18.9-20.4 

NAD 

828 

20.9-23.1 

Total  4,813 


*  Data  supplied  by  B.  Kyzer,  USAED,  Charleston. 
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Effects  of  Pinopolis  Inflow  on  Sediment  Sources  and  Shoallni 
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tt  Ocean  and  other  sediment  sources.  Assumed  to  be  related  to  estuarine  trap  efficiency.  See  paragraph  17  in  this 
appendix. 

f  ( 1  -  total  shoaling  volume/4.8)  *  100. 

^  Data  from  Patterson  (1983). 


APPENDIX  B:  FLOW  AND  SALINITY  MODEL 


I.  This  appendix  presents  and  describes  the  mathematical  equations  that 
make  up  the  flow  and  salinity  portions  of  the  FIBS  model  used  in  this  study. 
Due  to  the  complicated  nature  of  this  subject,  the  description  will  be  limited 
to  an  overview.  Blumberg  (1977)*  and  Wang  (1983  and  1984)  describe  similar 
laterally  averaged  flow  and  salinity  models.  Appendix  C  describes  the  sedi¬ 
ment  transport  model  equations,  Appendix  D  describes  the  fluid  mud  modeling, 
and  Appendices  E  and  F  contain  definitions  and  notation  used  throughout  the 
report  for  reference. 


Equations 

2.  The  following  governing  laterally  averaged  dynamic  differential  equa¬ 
tions  describe  estuarine  flow,  mixing,  and  circulation: 

a.  The  vertically  integrated  continuity  equation. 

b.  The  horizontal  momentum  equation. 

c.  The  continuity  equation. 

d.  The  salinity  transport  equation. 

Additional  terms  were  included  in  the  equations  to  describe  both  lateral  in¬ 
flows  and  tributaries.  The  following  algebraic  expressions  are  required  to 
dynamically  modify,  or  close,  the  set  of  differential  equations: 

a.  The  equation  of  state. 

b.  The  friction  coefficient  equation. 

£.  The  equations  for  vertical  eddy  dif fusivities  for  mass  and 
momentum. 

An  improved  vertical  eddy  diffusivity  closure  was  made  using  a  differential 
equation  for  the  transport  of  turbulent  kinetic  energy. 

3.  Laterally  averaged  equations  use  values  averaged  across  sections  at 
certain  increments  of  depth.  The  channel  widths  are  specified  at  each  node 
depth  as  indicated  in  the  hypothetical  channel  cross  section  (Figure  Bl). 


*  References  cited  in  this  appendix  can  be  found  at  the  end  of  the  main  body 
of  the  report. 
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Vertically  integrated  continuity  equation: 
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Horizontal  momentum  equation: 
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Continuity  equation  in  a  vertical  plane: 
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Salinity  transport  equation: 


ifri +  lx (SUB)  -  BKx 
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Equation  of  state; 


p  =  Pq  ( 1  +  Cl • S) 


where 

p  =  freshwater  reference  density  (0.9987  g/cu  cm) 

Cl  =  constant  (1.3751  E-6) 

S  =  salinity  parameter  expressed  here  as  chlorinity,  ppm 
Bed  friction: 
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2 

Cd(x)  *=  - SLT7?  (B6) 

2.22  R1'"5 

where  n  is  Manning's  coefficient  of  friction. 

Vertical  eddy  diffusivities: 

4.  Under  homogeneous  conditions,  characteristic  eddy  diffusivities  are 

N  =  K  =  K'  ,  where  K'  is  an  specified  constant.  The  effects  of  vertical 
zo  zo  z  z 

density  stratification  are  included  using  the  method  of  Munk  and  Anderson 
(1948).*  Thus: 

K  =  K  (1  +  3.33  Ri)"(3/2)  , 
z  zo 


and 


N 

z 


=  N 


zo 


(1  +  10  Ri)_(1/2) 
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*  References  cited  in  this  appendix  can  be  found  at  the  end  of  the  main  body 
of  the  report. 
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where 


l£  8P/3Z 

Po  ( 3U/ 0 Z) 2 


Then  at  depth. 


K  (z)  *  F  (z) 
z  z 


where 


is  a  similarity  distribution  for  vertical  eddy  diffusivity. 


Turbulent  kinetic  enert 


5.  Instead  of  assigning  a  value  to  ,  the  model  will  optionally  use 
turbulent  kinetic  energy  <  generated,  dissipated,  and  transported  by  the 
flow  to  calculate  vertical  diffusivities.  The  scheme  employed  is  similar  to 
that  described  by  Smith  and  Takhar  (1979)  and  Smith  (1982).  Transport  is 
computed  in  one  dimension  by: 


2 

3k  -  3k  ^  3  k 

75T  +  U  -s-v  -  D  — r 
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where  is  production  of  turbulent  kinetic  energy, 


pk  -  2ZA\r 


e  is  dissipation,  and 


e  -  1 


•»r) 


and  U  is  the  cross-sectional  average  velocity. 


0  -  i  f  MZ 


Then,  K  ■  N  ■  0.09  <  R  ,  and  the  values  for  K  ,  N  ,  K  (z)  and 
zo  zo  z  z  z 

N^(z)  are  computed  as  previously  described. 
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Boundary  Conditions 


Ocean  boundary 

6.  In  general,  the  varying  tidal  conditions  are  controlled  at  the  ocean 

boundary  of  the  model  mesh.  This  is  accomplished  by  specifying  a  tidal  water 

level  target  value  at  each  time-step  in  the  form  of  h'(t)  I  =  f(t)  .  The 

|  X—U 

model  has  the  capability  to  synthesize  tidal  sequences  for  any  specified 
starting  time  using  12  constituents  and  harmonic  coefficients  provided  by  the 
National  Ocean  Survey.  The  individual  constituents  used  included  M2  ,  S2  , 
K1  ,  01  ,  PI  ,  N2  ,  L2  ,  Sa  ,  M  f  ,  Mf  ,  Mm  ,  and  Ssa  . 

7.  A  nonreflecting  ocean  boundary  was  computed  for  this  study  by  re¬ 
placing  the  vertically  integrated  continuity  equation  with  the  following 
expression: 


where 


-^  +  Cw 
3t 


(B9) 
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x=o 


and  T^  is  a  damping  parameter  (found  to  be  about  500  sec).  This  method  is 
similar  to  that  presented  by  Blumberg  and  Kantha  (1985),  and  allows  long 
waves,  or  water  level  perturbations,  originating  within  the  model  domain  to 
pass  through  the  ocean  boundary  without  reflection.  The  non-reflecting  ocean 
boundary  permits  faster  stability  within  the  model,  and  in  this  study,  the 
model  reached  the  stable  condition  in  two  tidal  cycles.  During  the  flood 
tidal  phases,  salinity  concentration  was  specified  at  the  ocean  boundary; 
otherwise  the  transport  boundary  was  unconstrained. 

Upstream  boundary 

8.  Velocities  are  specified  at  upstream  inflow  boundaries  and  flow 
rates  are  specified  at  lateral  inflows.  In  this  study,  -he  condition 
U(t)  *  0  was  specified  at  Pinopolis,  the  upstream  end  The  model,  and  a 
lateral  inflow  was  specified  at  the  adjacent  element  to  represent  the  hydro- 
power  flow  release.  The  Cooper  River  contributes  almost  all  of  the  freshwater 
inflow  to  the  system  and  is  controlled  at  the  Pinopolis  Dam.  Salinity 
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concentrations  (usually  zero)  were  specified  at  these  inflows. 

Bottom  boundary 

9.  For  this  study,  the  bottom  boundary  (Z  =  0.0)  was  considered  to  be 
a  slip-flow  boundary.  The  shear  stress  was  imposed  by  the  equation: 

Nz  ||  -  Cd  U|uj  (BIO) 

Optionally,  the  model  will  extrapolate  shear  stress  linearly  from  within  the 
flow  to  the  bottom  boundary  using  a  method  suggested  by  Smith  (1982).  Within 
the  salinity  transport  equation,  no  flux  is  allowed  at  the  bottom  boundary 
and : 
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Vertically  Stretching  Coordinates 


10.  In  the  vertical  or  Z-direction,  computational  nodes  are  assigned  to 
the  surface,  intermediate,  and  bottom  levels.  When  changes  occur  in  the  water 
levels,  all  nodes  except  the  bottom  node  move,  or  stretch,  to  conform  with  the 
depth  of  flow.  The  typical  configuration  of  the  nodes  within  an  element  are 
shown  in  Figure  B2.  The  stretching  coordinates  system  introduces  new  terms 
into  the  governing  differential  equations.  A  derivative  in  some  level  of  the 
X-direction  is  calculated  by  the  following  expression: 
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The  X  and  Z  are  the  actual  world  coordinates.  The  x  and  z  are  the 
nondimensional  element  coordinates  that  vary  in  each  element  from  0.0  to  1.0 
upstream  and  vertically  upward,  respectively.  The  two  coordinate  systems  are 
related  by 
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Note  :  s/o/oes  are  c^uacdrat/c  (exaggerated  Aere). 

Figure  B2.  Example  computational  element  from  schematic  mesh 
showing  stretching  coordinates  and  nodes 
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e 

where  Ag  is  the  element  length.  This  method  is  similar  to  the  finite- 
difference  methods  presented  by  Sheng  (1983). 


(814) 


Solutions  of  the  Equations 


11.  The  governing  equations  were  solved  over  the  previously  described 
finite  elements,  using  an  explicit  orthogonal  collocation  method  of  weighted 
residuals.  Weighted  residual  methods  can  be  used  to  expand  an  unknown  solu¬ 
tion  in  a  set  of  trial  functions  T  over  a  domain  v  .  The  collocation 
method  applies  weighting  functions  (w^  ,  a  Dirac  delta  function)  to  the  trial 
functions,  chosen  to  be  set  of  orthogonal  polynomials.  The  node  locations 
Xj  satisfy  the  roots  of  the  polynomials.  The  form  of  the  collocation  method 
is  now 
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(  B 1 5 ) 

x . 

J 

where  w^  =  6  (x  -  ) . 

12.  Instead  of  solving  for  the  unknown  trial  function  coefficients,  a 
major  simplification  described  by  Finlayson  (1972)  is  that  the  solution  can  be 
derived  in  terms  of  values  at  the  collocation  nodes.  The  nodes  are  located  at 
Gauss  points  for  this  purpose,  and  the  method  is  known  as  orthogonal  colloca¬ 
tion.  The  derivatives  and  integrals  can  therefore  be  expressed  in  terms  of 
values  of  the  function  at  the  collocation  nodes  by 

dr  d2r 

g=Arand  ^  =  B  r  (B16) 

dx^ 

The  matrices  A  and  B  are  computed  from  the  trial  functions. 


Integration  in  Time 


13.  An  explicit  time  integration  method  was  used  to  compute  variables 
dynamically  over  time.  The  fourth-order  Runge-Kutta  method  was  used  in  this 
study.  This  method  is  very  accurate,  and  often  used  to  solve  "stiff"  equation 
sets,  equations  with  distinct  stability  limits.  Results  were  found  to  be 
independent  of  time-step  At  up  to  a  Courant  number  A^/AfCw  of  1.0  (or  2.0 
based  on  longitudinal  node  spacing).  With  the  initial  conditions  (t^.y^  over 
At  ,  the  function  y'(t)  =  f(t,y)  is  solved  in  four  steps: 

a. 

b. 

c. 

d. 


K„  =  At 


K,  =  At 
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Finally,  the  result  of  the  four  steps  are  combined  to  complete  the  time-step: 


yn+i  -  +  6  (Ko +  2Ki +  2k2  +  s) 


(B17) 
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APPENDIX  C:  SEDIMENT  TRANSPORT  MODEL  EQUATIONS 


1.  This  appendix  presents  and  describes  the  mathematical  equations  that 
make  up  the  sediment  transport  algorithms  for  fine-grain  sediment  within  the 
FIBS  model.  The  governing  differential  transport  equation  is  presented  first, 
with  the  remainder  of  this  appendix  briefly  describing  the  various  aspects  of 
sediment  transport  and  the  corresponding  mathematical  treatments. 

Characterization  of  Suspended  Sediment 

2.  Sediment  in  suspension  represent  some  specific  fraction  of  the  bed 
material  in  a  harbor  or  estuary.  Sediment  is  physically  sorted  or  fraction¬ 
ated  during  suspended  transport.  Fine-grained  sediments  are  hydraulically 
transported  almost  entirely  in  suspension  rather  than  as  bed  load.  Because  of 
the  differences  in  cohesion,  settling  characteristics,  etc.,  for  silts 

(4-72  pm)  and  clays  (0.45-4  pm),  f ined-grained  sediments  are  sometimes  charac¬ 
terized  as  a  sum  of  several  fractions  or  components.  The  fine-grained  mate¬ 
rial  can  also  contain  an  organic  fraction,  which  behaves  similar  to  cohesive 
sediments. 

3.  Fine-grained  sediments  exhibit  some  degree  of  cohesion;  thus  clay 
and  organic  solid  particles  aggregate  under  normal  estuarine  conditions.  The 
state  or  degree  of  aggregation  affects  the  deposition,  erosion,  and  settling 
processes,  described  later,  and  depends  on  sediment  concentration,  salinity, 
turbulence  in  the  flow,  pH,  temperature,  and  other  physical  factors. 

Governing  Transport  Equation 


4.  The  governing  differential  equation  for  laterally  averaged 
fine-grained  sediment  transport  is 
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where 


C  =  concentration  of  suspended  sediment 

W  =  settling  velocity 
8 

D  -  deposition 
E  =  erosion. 

Each  of  these  terms  is  described  by  a  set  of  algebraic  expressions,  and  is 
discussed  in  the  following  sections. 

Settling 

5.  Settling  is  that  component  of  suspended  particle  or  aggregate  motion 
caused  by  the  balance  between  gravity  and  viscous  drag  forces.  Settling  rates 
are  therefore  defined  in  quiescent  native  fluid.  Settling  characteristics 
affect  the  rates  of  deposition  and  the  vertical  distribution  of  suspended 
material  within  the  water  column. 

6.  Aggregation  is  very  important  to  cohesive  sediment  settling  rates, 
and  is  responsible  for  clay  deposition  in  estuaries  and  marine  environments. 
Aggregation  of  a  particular  sediment  particle  suspension  depends  primarily  on 
suspended  sediment  concentration,  current  shear  or  velocity  gradients,  and 
salinity.  Previous  experiments  on  the  effects  of  current  shear  on  settling 
found  impacts  at  shear  rates  above  those  encountered  within  most  natural  flows 
(Hunt  1982).*  Salinity  effects  on  aggregation  are  greatest  between  0-  and 
A-ppt  concentration.  Most  coastal  and  estuarine  harbors  are  almost  always 
above  this  range. 

7.  Three  ranges  of  concentration-dependent  settling  usually  occur.  At 
low  concentrations,  aggregate  and  particle  interaction  is  minimal,  and 
settling  is  independent  of  concentration.  At  intermediate  concentration, 
settling  is  enhanced  by  concentration  because  of  increased  aggregation  and 
particle  interaction.  At  high  concentrations,  aggregate  and  particle 
interaction  hinders  settling. 

8.  Suspended  sediment  concentrations  are  usually  in  the  enhanced- 

settling  range.  The  dependence  of  settling  velocity  W  in  the 

s 


*  References  cited  in  this  appendix  can  be  found  at  the  end  of  the  main  body 
of  the  report. 
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enhanced-settling  concentration  range  has  the  functional  form  (Ariathurai, 
MacArthur,  and  Krone  1977): 


W  =■  A1  CA2  (C2) 

s 

where  A1  is  a  constant  and  A2  is  the  enhanced  settling  exponent.  The  expo¬ 
nent  A2  is  usually  found  to  be  close  to  1.33.  The  concentration  range  over 
which  Equation  Cl  applies  varies  with  the  cohesive  properties  of  the  sediment. 
Generally  the  lower  bound  is  in  the  range  of  10-200  mg /£,  and  the  upper  bound 
is  in  the  range  of  2,000-75,000  mg /£. 

9.  Fine-grained  sediment  suspensions  usually  have  a  range  or  distri¬ 
bution  of  W  .  Clay  and  fine  silt  fractions  aggregate  to  form  a  relatively 
s 

unifora  settling  aggregate.  Medium  and  coarse  silt  fractions  settle  at  higher 
rates,  and  are  less  dependent  on  concentration  than  the  clay  fraction.  Set¬ 
tling  tests  can  be  used  to  determine  the  magnitude  and  distribution  of  at 

various  suspended  sediment  concentrations  for  the  finer  fractions  of  the 
material.  Such  testing  is  done  in  the  field  or  laboratory,  or  a  combination 
of  both. 


Deposition 

10.  Deposition  D  ,  or  flux  of  sediment  material  to  the  bed,  is  the  sum 
over  a  number  of  fractions  of  settling  flux  times  deposition  probability 
(Mehta  et  al.  1986) 


k 

D  =  V*  p  -W  C.  (C3) 

Z-#  7  s  J 

j-1  J 


where 

k  *  number  of  sediment  fractions 

P  •  probability  that  an  aggregate  which  has  reached  the  bed  will  remain 
remain  there 

C  ■  concentration  just  above  the  bed 
j  -  subscript  indicating  a  sediment  fraction 
P  varies  linearly  from  0  at  a  bed  shear  stress  equal  to  the  critical  shear 
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for  deposition,  “  tc^  ,  to  1  at  zero  bed  shear  stress,  =  0  .  The 
functional  form  1  -  (t^/t^)  where  <  tc<j  is  used  for  P  (Krone  1962). 
Laboratory  deposition  testing  can  be  used  to  determine  rccJ  and  the  magnitude 
of  the  product  PWg  for  each  sediment  fraction  identified. 

11.  A  suspension  of  uniform  material  in  a  steady,  uniform  flow  will 
either  deposit  completely  or  remain  entirely  suspended  depending  on  whether 

T,  is  below  or  above  ,  according  to  Equation  C3.  The  consequence  of  the 

presence  of  multiple  sediment  fractions  in  a  suspension  is  that,  under  a 
given  flow  condition,  some  sediment  fractions  may  deposit  while  others  may 
remain  in  suspension.  The  suspension  may  therefore  transport  an  equilibrium 
concentration  (some  fraction  of  the  source  concentration)  indefinitely. 

12.  The  values  of  W  inferred  from  deposition  tests  are  smaller  than 

s 

those  obtained  from  quiescent  settling  tube  tests.  The  cause  for  this  is  not 
known.  However,  shear  in  the  flow  is  greatest  just  above  the  bed,  and  could 
cause  disaggregation  and/or  produce  lift  forces  counteracting  settling  at  this 
point. 


Re  suspension 


13.  Resuspension  of  bed  material  occurs  in  four  modes: 

a.  Redispersion 

b.  Particle  erosion 

c.  Significant  erosion 

d.  Mass  erosion 

Resuspension  is  the  most  difficult  fine-grained  sedimentation  process  to 
either  test  or  form  predictions  for  occurrence.  Redispersion  occurs  when  the 
bed  is  fluid  and  particles  are  fluid  supported,  lacking  interparticle  or 
interaggregate  structure.  Redispersion  proceeds  as  an  entrainment  process, 
scaled  to  density  difference  between  the  fluid  mud  and  overlying  water,  an 
appropriate  velocity  or  shear  velocity,  layer  thickness,  and  viscous  effects. 
Redispersion  is  initiated  at  low  bed  shear  stresses  and  is  generally  rapid. 
Mass  erosion,  on  the  other  hand,  has  been  postulated  to  occur  at  very  high 
,  and  is  initiated  when  a  subbottom  layer  fails  under  extreme  hydraulic 
shear  conditions,  allowing  a  large  bed  layer  to  be  carried  away  by  the  flow. 
Mass  erosion,  under  these  conditions,  is  not  well  documented,  and  not  expected 
to  occur  under  normal  or  even  severe  estuarine  conditions.  Mass  erosion  can 
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also  occur  in  quasi-uniform  fluid  mud  layers  on  a  slope,  where  stresses  in¬ 
crease  vertically  downward  and  exceed  the  shear  strength  of  the  layer,  and  the 
layer  slides  as  a  mass. 

14.  Particle  and  significant  erosion  modes  are  similar.  At  t,  above 

b 

a  critical  value,  particles  or  clusters  of  particles  are  individually  dis¬ 
lodged  from  the  sediment  bed  as  interaggregate  bonds  are  broken.  Particle 
resuspension  is  related  to  the  shear  stress  in  excess  of  a  critical  value,  and 
to  an  erosion  rate  constant  M  (Ariathurai,  MacArthur,  and  Krone  1977): 


E  *  M 


(C4) 


where 

E  =  erosion  rate 

=  critical  erosion  shear  stress  for  particle  erosion 
Observed  erosion  does  not  follow  Equation  C4  indefinitely.  Suspension  con¬ 
centrations  above  experimental  eroding  beds  often  reach  constant  or  equilib¬ 
rium  values  which  depend  on  a  combination  of  the  bed  shear  stress  and  char¬ 
acter  of  the  bed.  Equilibrium  fine-grained  suspensions  form  as  erosion  rates 
decrease  with  time  to  zero,  while  the  flow  remains  constant.  Equilibrium 
fine-grained  suspensions  have  been  found  to  be  related  not  to  the  transport 
capacity  of  the  flow  (as  for  sand),  but  to  vertical  differences  or  nonhomo¬ 
geneity  in  the  bed  (either  particle  characteristics  or  bed  density)  or  to 
armoring  by  selective  erosion  at  the  bed  surface.  Resuspension  tests  can  be 
used  to  determine  the  magnitude  of  M  and  for  representative  sediment 

fraction,  and  to  detect  the  formation  and  nature  of  equilibrium  suspensions. 

15.  Typical  laboratory  relationships  between  E  and  t  show  two 

b 

regimes  or  modes  of  erosion  (Hunt  1981).  At  low  shear  stresses,  erosion  has 
been  found  to  proceed  first  as  surface  particle  flaking  which  has  a  low 
threshold  and  a  lower  value  of  M  .  At  higher  values  of  ,  a  break  point 

is  passed  at  which  E  and  M  increases  sharply.  This  is  termed  significant 
erosion,  and  was  first  observed  by  Partheniades  (1962).  Visual  observations 
by  the  author  have  suggested  that  significant  erosion  proceeds  as  small  bed 
sediment  chunks  are  detached  by  the  smallest  scale  turbulence  (with  character¬ 
istic  lengths  on  the  order  of  a  millimetre).  At  these  higher  levels  of  t,  , 

b 


C5 


the  smallest  eddies  have  diminished  size  and  increased  energy.  During  turbu¬ 
lent  bursting  of  the  laminar  sublayer,  small  eddies  impinge  on  the  bed  and 
dislodge  chunks  of  about  their  own  size.  Extrapolating  the  higher  erosion 
mode  yields  a  relatively  high  apparent  critical  erosion  threshold.  Figure  Cl 
shows  an  example  plot  of  the  general  form  described. 

16.  Significant  erosion  is  described  similarly: 


E 


Tch 


T 

C 


(C5) 


where 


“h 

Tch 


the  significant  erosion  rate  constant 

the  critical  shear  stress  for  significant  erosion 
(t  k  Is  larger  than  ,  and  is  larger  than  M) 


-7J  ,  -V/8?  m 


Figure  Cl.  Typical  nonlinear  erosion  behavior  of 
cohesive  sediment  (data  from  Hunt  1981) 
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APPENDIX  D:  FLUID  MUD  MODELING 


Background 


1.  The  Fine-grained  Bed  pediment  or  FIBS  model  was  formulated  to 
address  estuarine  conditions  where  beds  are  soft  muds.  Fine-sediment  bed  con¬ 
ditions  are  not  uniquely  related  to  grain  size,  but  are  time-dependent,  and 
related  to  admixture  composition  and  characteristics.  Therefore,  bed  condi¬ 
tions  depend  on  conditions  within  the  model  domain,  and  are  properly  an  inter¬ 
nal  rather  than  external  model  quantity.  The  sections  that  follow  describe  a 
dynamic  fluid  mud  model. 

2.  Fluid  mud  has  concentrations  above  that  at  which  the  maximum  set¬ 

tling  flux  occurs,  and  below  that  at  which  important  sedimentary  structure 
forms.  Settling  flux  is  the  product  of  settling  velocity  and  concentra¬ 

tion.  Settling  velocities  decrease  with  increased  concentration,  and  a  con¬ 
centration  is  reached  where  settling  flux  (the  product  of  concentration  and 
hindered  settling  velocity  W^)  begins  to  decrease.  This  point  is  the  maximum 
settling  flux,  and  occurs  generally  at  concentrations  between  5  and  100  g/1 
for  estuarine  cohesive  sediments.  The  behavior  of  suspensions  changes  rather 
abruptly  at  this  concentration,  and  suspensions  begin  to  settle  as  a  mass.  A 
clear  layer  forms  above  a  distinct  interface.  The  descent  of  that  interface 
is  associated  with  consolidation,  which  is  used  here  to  denote  volume  change 
by  the  effect  of  settling. 

3.  The  settling  characteristics  of  dense  suspensions  can  be  determined 
by  column  tests  conducted  over  a  range  of  concentrations.  Results  are  fit  to 
any  number  of  functional  relationships  and  used  in  further  numerical  analysis. 
In  fluid  mud,  settling  of  sediment  particles  is  hindered  by  particle  interac¬ 
tion  and  by  interaction  between  particles  and  the  pore  water.  The  higher  the 
fluid  mud  concentration,  the  more  slowly  the  mixture  consolidates  by  settling. 
Settling  velocities  decrease,  sometimes  sharply,  with  increasing  initial 
suspension  concentration. 

4.  Several  functions  can  be  used  to  fit  high-concentration  settling 
data.  The  best  known  functional  descriptor  of  the  effect  of  concentration  on 
hindered-settling  velocity  appears  to  be  the  Richardson-Zaki  (1954)*  equation: 

*  References  cited  in  this  appendix  can  be  found  at  the  end  of  the  main  body 
of  the  report. 
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Wh  -  W1  ( 1-B2  Cg)B1  (Dl) 

where  is  the  hindered-settling  velocity,  W1  is  a  reference  settling 

velocity,  B2  is  the  inverse  of  the  fully-settled  concentration,  and  C  is 

s 

the  dense  suspension  concentration.  Richardson  and  Zaki  found  the  exponent  B1 
to  have  a  value  of  4.7.  Sediments  from  several  United  States  estuaries  have 
been  found  by  the  author  to  fit  the  Richardson-Zaki  equation  reasonably  well. 
However,  not  all  data  do  so  (Owen  1970). 

Yield  stress 

5.  Rheological  tests  suggest  that  fluid  mud  is  a  viscoelastic  material 
with  the  ability  to  sustain  a  finite  shear  stress  with  zero  continuous  strain. 
The  yield  stress  for  a  number  of  estuarine  muds  was  determined  by  Krone 
(1963),  who  assumed  a  Bingham  model  and  found  yield  stress  to  be  strongly 
dependent  on  concentration.  Rheological  measurements  on  20-  to  60-g/l  suspen¬ 
sions  with  a  rotating  cylinder  viscometer  were  found  to  obey  a  power  law  with 
a  slope  of  2.5.  The  yield  stress  of  a  sediment  material  is  directly  related 

to  the  critical  shear  stress  for  erosion  T  ,  so  that 

c 

Tc  -  Dl  CsD2  (D2) 

where  C  is  the  concentration  of  fluid  mud  as  dry  weight  of  sediment  per 
s 

unit  volume  of  mud,  and  Dl  and  D2  are  imperical  constants.  Thus,  the  hydrau¬ 
lic  shear  strengths  of  dense  suspensions  with  24  g/1  and  60  g/1  solids  varied 
by  about  a  factor  of  10. 

6.  Likewise,  the  significant  erosion  yield  or  critical  stress  may  be 
strongly  dependent  on  bed  concentration. 

T  .  *  Dhl  C  Dh2  (D3) 

cn  s 

where  Dhl  and  Dh2  are  empirical  constants  for  significant  erosion. 

7.  The  erodibility  of  estuarine  sediments  has  also  been  observed  to 
depend  on  concentration  as  indicated  by  Equation  D3.  Teeter  (1987)  found  the 
exponent  Dh2  to  be  3.38  for  San  Francisco  Bay  sediments.  The  experimental 
concentration  range  for  the  bed  was  from  275  to  525  g/1. 
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Model  Description 


8.  A  transport  theory  for  the  consolidation  by  hindered  settling  can  be 
constructed  using  a  form  of  the  vertical  transport  equation  with  diffusion  and 
by  nondimensionalizing  depth  as  z  =  Z/Hg  : 
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where  H  is  the  total  depth  or  height  of  the  dense  suspension,  and  K  is  a 
s  s 

diffusivity  for  mass.  Diffusive  transport  in  fluid  muds  could  come  from  bio- 
turbation  (on  the  order  of  IE-7  sq  cm/sec,  for  example),  or  from  convective 
instabilities  observed  in  actively  settling  suspensions.  Numerically,  the 
diffusion  term  was  employed  to  counteract  instabilities  induced  by  the 
boundary  conditions  required  for  high  deposition  rates. 

9.  Kynch  (1952)  advanced  a  kinematic  theory  for  the  analysis  of 
hindered-settling  tests  based  on  a  continuity  equation  similar  to  Equation  DA. 
He  assumed  that  settling  velocity  was  a  local  characteristic  dependent  on  the 
concentration  of  the  suspension.  In  the  development  of  his  theory,  Kynch 
pointed  out  that  the  total  derivative  of  concentration  (he  used  density)  for  a 
moving  reference  point  in  a  dense  suspension  is  equal  to  the  partial  deriva¬ 
tive  of  concentration  with  respect  to  time  plus  the  convection  of  the  concen¬ 
tration  gradient.  Thus  in  the  present  notation: 
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10.  Discontinuities  or  step  changes  in  the  concentration  field  compli¬ 
cate  the  analysis  of  dense  suspensions.  The  interface  observed  during  consol¬ 
idation  experiments  is  such  a  discontinuity.  Kynch  (1952)  explained  how 
several  transport  conditions  can  produce  discontinuities.  Discontinuities  are 
associated  with  the  concentration  at  which  the  vertical  flux  is  maximum,  and 
other  points  where  the  rate  of  change  in  vertical  flux  with  respect  to  concen¬ 
tration  is  zero.  The  former  condition  always  occurs,  while  the  latter  can 
occur  if  the  relationship  between  concentration  and  settling  velocity  is  not 
smooth. 

D3 


11.  Solution  of  the  model  governing  equations  requires  that  the  differ¬ 
ential  equations  be  replaced  by  matching  conditions  at  concentration  discon¬ 
tinuities.  The  location  of  discontinuities  therefore  becomes  part  of  the 
problem  and  must  be  solved  by  the  model.  Numerical  solution  of  transport 
equations  apparently  cannot  span  concentration  discontinuities  without  immedi¬ 
ately  becoming  unstable.  Numerical  methods  used  to  solve  hindered  settling 
'-^tinuity  equations  were  very  similar  to  those  used  for  flow  and  salinity. 

12.  The  model  mesh  was  contracted  during  computations  as  the  dense  sus¬ 
pension  collapsed  so  that  z  =  1.0  was  fixed  to  the  upper  boundary.  The 
height  of  the  suspension  to  the  interface  or  discontinuity  was  calculated,  for 
instance,  at  any  time  by  conservation  of  mass  as 


dH 

— —  =  -W 
dt  wh 


(D6) 


where  the  subscript  i  indicated  the  location  of  the  interface.  A  no-flux 
bottom  boundary  condition  W^C  =  0  was  imposed  at  z  =  0  .  No  condition  was 
required  at  the  upper  boundary  of  the  dense  suspension  since  the  equation  was 
applied  only  from  fixed  bed  to  the  moving  interface  (z  =  1.0).  Initial  height 
and  concentration  of  the  dense  suspension  were  starting  points  for  model 
computations . 

Depositional  formation  of  dense  suspensions 

13.  The  concentration  of  newly  deposited  material  affects  the  thickness 
and  hydraulic  shear  strength  of  the  deposit.  The  concentrations  of  material 
as  they  undergo  deposition  from  suspension  and  form  a  sediment  bed  have  been 
studied  very  little  (Teeter  1986). 

14.  During  depositional  formation,  depositional  flux  C  Wg  must  be 
included  in  the  model  dense-suspension  height  equation: 


dH  C  W 

dT  "  -c-*  ’  wh  (D7) 

si  i 

where  C  is  the  suspension  concentration.  However,  at  high  deposition  rates 
this  equation  does  not  reproduce  the  observed  fluid  mud  formation  very  well. 
Figure  Dl  shows  a  plot  of  the  experimentally  determined  variables  which  appear 
in  the  equation,  and  shows  a  disagreement  between  the  observed  fluid  mud 
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Figure  Dl.  Example  experimental  data  on  the  terms  of 

Equation  D7 

characteristics  and  the  equation.  Rapid  deposition  apparently  reduces  hin¬ 
dered  settling  in  the  deposit.  At  higher  concentrations  and  settling  rates 

the  term  C  W  /C  was  closer  to  (dH  /dt)  -  W,  .  Therefore,  a  criterion 
s  s_^  s 

was  used  to  define  the  higher  depositional  range,  and  that  criterion  was  found 

to  be  dH  /dt  >  2  W, 
s  h . 

1 

15.  The  upper  boundary  condition  imposed  on  the  transport  equation  was 
changed  during  rapid  deposition  to 


C  =  B3  W  (D8) 

Si  Si 


where  113  is  an  empirical  constant.  Figure  D2  shows  a  plot  of  Wg  versus 
initial  for  tests  performed  with  Corpus  Christi  sediments. 

Verification  to  laboratory  experiments 

16.  The  hindered-settling  consolidation  model  has  reproduced  observed 
interface  descent  for  several  sediments.  Laboratory  data  from  hindered 
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Initial  Suspension  Concentration 


high 


medium 


W3  ,  mm/see 

A  C;  determined  at  peak 
deposit  thickness 

□  C7  determined  before  peak 
deposit  thickness 

Figure  D2.  Example  experimental  data  on  the 
relationship  between  C  at  the  interface 
of  the  dense  suspension,  and  W  of  the 
overlying  suspension 

settling  consolidation  tests  were  used  to  fit  Equation  Dl.  Figure  D3  shows  an 
example  comparison  between  numerical  and  experimental  data  for  Atchafalaya  Bay 
sediments.  The  computed  vertical  concentration  profiles  were  also  very  simi¬ 
lar  to  those  observed  during  tests.  Concentrations  increased  at  the  bottom 
first.  The  commonly  observed  inflection  in  the  interface  descent  occurred  in 
the  numerical  simulations  when  the  surface  concentration  first  increased.  The 
rate  of  interface  descent  was  found  not  to  be  sensitive  to  the  initial  height 
of  the  suspension,  as  reported  by  several  investigators. 

17.  Results  of  numerical  calculations  of  fluid  mud  formation  are  shown 
in  Figure  D4.  Both  the  dense  suspension  and  overlying  suspension  were 
modelled  using  the  previously  described  method  for  rapid  deposition.  Note 
that  deposition  ceased  at  about  2.2  hours  into  the  test  as  the  overlying  sus¬ 
pension  became  exhausted  of  8*  ment.  The  dense  suspension  deposit,  which  had 
grown  linearly  with  time,  began  to  consolidate  at  this  point. 
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Figure  D3.  Example  comparison  between  experimental  data  (■)  and  model 
results  for  the  interface  height  of  a  hindered-settling  consolidation 

test 
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APPENDIX  E:  DEFINITIONS  OF  TERMS 


Bed 

Deposition 

Flux 

Flux  components 

Resuspension 

Salinity 

Settling 
Stokes  velocity 

Stratification 

Tidal  average 

Tidal  hydraulics 

Vertical  circulation 


The  lower  level  of  the  flow  field,  i.e.,  the  level 
of  no  flow 

Removal  to  the  bed  of  sediment  from  suspension  by 
settling. 

The  transport  of  salt  or  suspended  sediment  through 
a  certain  area.  Calculated  as  the  product  of  veloc¬ 
ities,  salinities,  or  suspended  sediment  concentra¬ 
tions,  and  cross-sectional  or  unit  areas,  and 
usually  summed  over  a  tidal  cycle. 

Statistical  correlations  calculated  by  a  sequence  of 
time  and  depth  averaging  for  the  purpose  of  resolv¬ 
ing  instantaneous  and  depth  deviations  in  fluxes  at 
a  sampling  station. 

Erosion.  The  removal  of  sediment  from  the  bed  by 
the  forces  of  the  flow. 

Concentration  by  weight  (expressed  as  parts  per 
thousand  or  ppt)  of  inorganic  matter  (mainly  chlo¬ 
ride,  bromide,  sodium,  potassium,  magnesium,  and 
calcium)  in  seawater  or  brackish  water  (dilute 
seawater) . 

The  balance  between  gravitational  and  viscous  forces 
on  suspended  particles. 

A  residual  flow  generated  by  tide  or  other  long-wave 
propagation.  The  difference  between  Eulerian  and 
Lagrangian  tidal-averaged  velocities. 

Vertical  salinity  and  density  distributions  that 
stabilize  estuarine  flows  by  buoyancy  effects,  and 
that  inhibit  vertical  turbulent  transport  of  salin¬ 
ity,  suspended  matter,  and  momentum. 

The  average  of  a  quantity  or  function  over  a  tidal 
cycle.  In  the  case  of  flows  or  fluxes,  the  residual 
of  tidal  motion. 

Instantaneous  flows  and  water-surface  elevations 
associated  with  earth/astronomical  gravitational 
effects . 

Tidal  residual  flow  in  the  vertical  plan  charac¬ 
terized  in  estuaries  by  upstream  flow  at  the  bottom 
and  seaward  flow  near  the  surface. 


El 


Vertical  mixing  The  vertical  turbulent  exchange  of  salinity  and 

suspended  material  that  is  generated  by  an  estuarine 
flow. 
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APPENDIX  F:  NOTATION 


Subscripts  and  indices: 
b  Bottom;  bed 

f  Fresh  water 

i  Quantities  at  the  suspension/bed  interface 

j  Sediment  fraction  index 

l  Lateral  inflow  quantaties  such  as  salinity 
o  Indicates  a  depth  average  or  vertically  homogeneous 
t  Time 

v  Vertical  deviation 

x  Horizontal  coordinate 

z  Depth  or  vertical  coordinate 

Parameters: 

B  Breadth  (width) 

Bq  Breadth  (width)  at  water  surface 

C  Concentration  of  suspended  sediment 

Cd  Quadratic  friction  coefficient 

C  Concentration  of  bed  sediments 

s 

Cw  Shallow-water  wave  speed 

D  Deposition 

One-dimensional  unadjusted  diffusion  coefficient 
E  Erosion 

H  Instantaneous  water  depth 

Hg  Thickness  of  a  fluid  mud  or  dense  suspension 

K  Horizontal  and  vertical  eddy  diffusivity  for  mass  coefficients 

X  f  z 

M  Erosion  rate  constant  for  particle  erosion 


FI 


Erosion  rate  constant  for  significant  erosion 

N  Horizontal  and  vertical  eddy  diffusivity  for  momentum  coefficients 
X  f  z 

P  Probability  of  erosion 

P^  Production  of  turbulent  kinetic  energy 

R  Water  depth  to  mean  tide  level 
Ri  Richardson  number 
S  Salinity 

Salinity  of  lateral  inflow  or  inflowing  from  a  branch 

U  Horizontal  (laterally  averaged)  velocity 

W  Vertical  (laterally  averaged)  velocity 

Hindered  settling  velocity 

W  Settling  velocity 
8 

g  Acceleration  due  to  gravity 

h  Deviation  of  the  water  surface  from  the  mean  tide  level 
k  An  index,  such  as  the  number  of  sediment  fractions 
n  Manning's  friction  coefficient 
q  Lateral  inflow  rate 

x  Nondimensional  horizontal  coordinate  for  X 

z  Nondimensional  vertical  coordinate  for  Z 

ic  Turbulent  kinetic  energy 
p  Density 

Po  Reference  density  (freshwater) 

A  Element  length 

e 

t  Shear  stress 

e  Dissipation  of  turbulent  kinetic  energy 
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